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Abstract 

We review the status of the practical operator product expansion (OPE), when apphed to two- 
point correlators of QCD currents which interpolate to mesonic resonances, in view of the violations 
of local quark-hadron duality. Covered topics are: a mini-review of mesonic QCD sum rules in 
vacuum, at finite temperature, or at finite baryon density, a comparison of model calculations 
of current-current correlation functions in 2D and 4D with the OPE expressions, a discussion of 
meson distribution amplitudes in the light of nonperturbatively nonlocal modifications of the OPE, 
and a reorganization of the OPE which (partially) resums powers of covariant derivatives. 
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1 Preface 



We start by giving a short outline of the present review article: 

In the next five sections we lay the foundations for reviewing (some of) the research on local quark- 
hadron duality violations being conducted within the last twenty years or so. We introduce the operator 
product expansion (OPE) as an expansion of a current- current correlator in powers of inverse Euclidean, 
external momenta and the strong coupling constant as- We distinguish between the full OPE, which 
most probably is an asymptotic series, and truncations thereof used in practice. A definition of local 
quark-hadron duality, which conceptually rests on the OPE, is given. The concept of the OPE is 
applied to hadronic physics in the framework of the Shifman-Vainshtein-Zakharov (SVZ) sum rules 
associated with light- and heavy-quark mesons in vacuum and in a hot or baryon rich medium. We do 
not discuss the sum-rule program for baryonic resonances. We also review briefiy the relation between 
power corrections and renormalons. Sees. 2 to 5 naturally have large overlaps with existing review 
articles and books on QCD sum rules and on theoretical as well as on practical aspects of the OPE 

EDinni. 

In Sec. 6 we take a closer look at violations of local duality. We first compare the predictions of 
analytically continued practical OPEs with the experimental data. Second, we review an instanton 
model calculation of the light-quark current-current correlator in Euclidean spacetime and compare the 
result with the practical OPE of this correlator. Third, we review an analysis of the current-current 
correlator in the 't Hooft model. The focus is on oscillatory components in the decay width of heavy 
mesons (as a function of the meson mass) when taking 0{1/Nc) corrections into account. The adopted 
point of view in the latter two investigations is that the model calculations are in some sense realistic, 
that is, they largely resemble the experimental situation. Disagreement between the prediction of 
hadronic spectra resting on the model calculation on the one hand and the OPE on the other hand 
thus are interpreted as violations of local quark-hadron duality. In Sec. 16. II we give an overview on 
experimental hadron spectra induced by the electromagnetic current and by r decays. We also address 
Bg-Bg mixing: a problem which can be tackled by appealing to local quark hadron duality. A rather 
large overlap of Sec. 16. 21 and Sec. 16. 31 with the reviews [32j exists. In the present review the analysis in 
the 't Hooft model is presented in a more self-contained way than in |32j . 

A discussion of the effects of nonperturbative nonlocality in modified practical OPEs is carried out 
in section 7 of the review. In a first step the relation between power corrections in an OPE and mesonic 
distribution amplitudes is discussed. Subsequently, the theoretical and phenomenological need for the 
inclusion of nonlocal condensates is pointed out. We also review the systematic inclusion of nonpertur- 
batively nonlocal quantities in the OPE-based description of the hadron-to-hadron, hadron-to- vacuum, 
and vacuum-to- vacuum matrix elements of time-ordered current-current products. Phenomenologically 
consequences are explored in the latter case. To the best of the author's knowledge no review article 
exists which would have a substantial overlap with section 7. 

The presentation is not very technical. Rather, we have tried to capture the basic theoretical 
concepts and their applications. Most of the time our discussion aims at a direct comparison with the 
experimental situation. This review is by no means complete - our apologies to all those authors whose 
contributions are not explicitly referred to. 

2 Introduction 

Quantum Chromodynamics (QCD), the non-Abelian gauge theory of interacting quarks and gluons, is 
by now the widely accepted microscopic theory of strong interactions. Embedded into the Standard 
Model of particle physics it has passed various experimental tests. The expansion about the situation of 
asympotically free fundamental degrees of freedom at large momenta or particle masses is a conceptually 
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and practically appealing feature P3I21IS1- It has a wealth of applications: the perturbative calculability 
of the evolution of hadron structure functions [H El El 13' ^ perturbative description of quark or gluon 
induced jets in e^e~ annihilations jH], the perturbative matching of the full theory to effective theories 
for heavy flavors El El El C31 El El CSI and the matching of the electroweak sector of the Standard 
Model involving perturbative QCD corrections with a low-energy effective theory describing the weak 
decays and mixings of heavy flavors, see for example ^3 Ej for 6 — ^ s decays. Last but not least, 
at large external momenta asymptotic freedom guarantees the usefulness of an expansion of Euclidean 
correlation functions of hadron- interpolating currents starting from the parton model fniEDj- This is 
crucial for the method of QCD sum rules [2111221 whose explanation is the starting point of the present 
review. 

While QCD is well understood and tested at large external momenta p ^ 0.5 GeV the strongly 
coupled low-energy regime denies a perturbative treatment. The only presently known fundamental 
approach to learn about the dynamics of the relevant degrees of freedom populating the ground state 
of QCD at low energy and being responsible for the hadronic spectrum and hadron dynamics are 
lattice simulations [201 • Analytical approaches are not clear-cut - not even the identification of the 
relevant degrees of freedom is unique -, and explanations of low-energy phenomena have only been 
found partially so far. For example, the instanton j2S| gas provides a plausible, microscopic mechanism 
for dynamical chiral symmetry breaking [23j and the large mass of the rj' |21| but unless interactions 
between instantons are parametrized as in the instanton-liquid model (see [211 foi' a review) it does not 
explain color confinement. Dual models j2Zl based on the idea that confinement is realized through 
the dual Meissner effect |28| 1^ explain confinement almost by definition and can not be considered 
fundamental. 

QCD sum rules provide a successful analytic way to 'scratch' into the nonperturbative regime. Many 
exhaustive reviews of the field exist [HI] , and therefore we will limit ourselves to a very basic introduction 
in the next section. The sum-rule approach is pragmatic in the sense that it does not aim at calculating 
nonperturbative ground-state parameters or hadronic wave functions from first principles. Rather, 
QCD sum rules extract universal nonperturbative parameters from limited experimental information to 
subsequently use them for predictions. The sum-rule method has been extended to investigate hadronic 
properties in the presence of a hot [73| or a dense hadronic medium (for a review see |lU8j ) which is of 
great relevance for ongoing experiments with colliding (ultra)relativistic heavy ions. 

The basic object of investigation in QCD sum rules is an A^-point correlator of QCD currents. The 
sum-rule method rests on analyticity in the external momentum variable (s) and on the assumption of 
quark-hadron duality. In its strong, local form, quark-hadron duality is the situation that a hadronic 
cross section can be related pointwise to a theoretical expression obtained in terms of quark and gluon 
variables. At low energy strong limitations on the accuracy of the theoretical expression exist. One 
possible theoretical approach is the analytical continuation of the OPE 36 which is an asymptotic 
series in powers of the strong coupling and in powers of the inverse external momentum (up to 
logarithms) [211 1142j . To circumvent the problem associated with the asymptotic series one appeals to 
analyticity in the external momentum variable to relate the incompletely known theoretical part to an 
average over the hadronic cross section by means of a dispersion relation [21] . This pragmatic approach 
is rather fruitful. 

Different meanings are attributed to quark-hadron duality in the literature. For the present review, 
which is mainly concerned with its local form, we use the definition given by Shifman in |144j for a 
one-variable situation such as the e"'"e~ annihilation into hadrons. Shifman's definition is as follows: 
The current correlator, associated with the process of hadron creation out of the vacuum, is theoret- 
ically evaluated in terms of quark and gluon fields at an external momentum q with = —Q^ < 
(Euclidean region) using an OPE Approximating this current correlator by a truncation of the 
purely perturbative part of the OPE at and a truncation at power (not counting logarithms 
arising from anomalous dimensions of contributing operators), the uncertainty of the result should be of 
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orders a^J'^^^ and Q^^^^^^ if the OPE represents a converging series. In practice one usually has k = 2 
or k = 3 and D = 6 or D = 8 (practical OPE). The uncertainty in the truncated OPE would translate 
into an uncertainty of order 0;^'^+^) and s~(^+i)/2 of the theoretically predicted spectral function Pthool-s)- 
The latter is obtained by calculating the imaginary part of a term-by-term analytical continuation to 
negative Q^, —Q^ = s > 0, of the truncated OPE. If the experimentally measured spectral function 
Pexp(s) coincides with pthco(s) up to the above uncertainty within a certain range of s values we say that 
the quark-gluon prediction, resting on the OPE, is dual to the hadronic spectral function in this range. 
If this is not true we speak of a violation of local quark-hadron duality. While local quark-hadron 
duality is badly violated in practical OPEs at momenta close to the lowest resonance mass there usually 
is a window of Euclidean momenta Q where the practical OPE is seen to be equal (or dual) to the 
dispersion integral over the spectral function. We will refer to this kind of duality as global or sum-rule 
duality. It is known for a long time that in some exceptional channels QCD sum rules do not work 
[SHI ISH ISni HO] ■ Namely, a stability under the variation of quantities parametrizing our ignorance about 
the hadron spectrum and the vacuum structure does not exist in these channels. 

3 SVZ sum rules 

To prepare for the discussion of duality violation in later chapters we briefly review the method of 
QCD sum rules of Shifman, Vainshtein, and Zakharov (SVZ)|2T1 122] focussing on mesonic two-point 
correlators. A remark on the notation is in order: In this and subsequent sections multiple meanings 
are attributed to the symbol T. It may stand for a correlator or time-ordering or the temperature. The 
symbol T is the inverse Borel transform of B[T]. What is meant should become clear from the context. 
We apologize for the inconvenience. 

Let us start by giving a more specific characterization of the SVZ approach than in the last section: 
The idea is to express hadron parameters in terms of expectation values of gauge invariant operators 
which are composed of fundamental quark and/or gluon fields. In practice, the hadron parameters 
of the lowest states are studied. The lowest states are explicitly presented in models for the spectral 
functions. The higher resonances are absorbed into a continuum which is usually approximated by 
perturbative QCD. To separate the continuum from the lowest state an effective threshold sq, which is 
a model parameter, is introduced. The parameter sq is to be determined by the sum rule by requiring 
hadron parameters to be least sensitive to its variation. Obviously, one does not assume local quark- 
hadron duality here - only the integral over the perturbative QCD spectrum is assumed to be equal to 
the integral over the experimental spectrum down to the threshold Sq. 

The basic object of investigation is the vacuum correlator of two gauge invariant, hadron-interpolating 
currents jl and j^, in momentum space where L, L' and F, F' denote Lorentz and flavor quantum num- 
bers, respectively 

T[[:{q) ^zjd'x e^^^' (o|T {j[{x)jI: (0)) |o) . (1) 

Quantum numbers with respect to discrete symmetries are suppressed. The quantum numbers of a 
currents correspond to the probed hadronic resonance or the probed mixing of them. In Eq. (^Q) T refers 
to the time-ordering symbol. Notice that the mass-dimension of the correlator on the left-hand side of 
Eq. (0) is two if the currents are quark bilinears such as 

j; = V'y7MV^, ^ = {u,df, (z = l,...,3). (2) 

Here r* stand for SU(2) generators in the fundamental representation normalized to TrrV-' = 2 6^^ . If 
a vector current j'J is conserved then a current conserving (transverse) tensor structure can be factored 
out of its correlator, and we have 

T^f{q) = {q,q.-q'9,u)T^^{q'). (3) 



5 



■ • 



-7 Req^ 



Figure 1: The integration contour in the complex plane. The dashed line indicates the 
part at infinity. Dots denote isolated poles related to stable particles. The thick line is the 
continuous part of the spectrum. 



In the case of a quark-bilinear current this implies the scalar amplitude T^^(g^) to be a dimensionless 
function of the square of the external momentum q^. In the more general case of non or partially 
conserved currents, such as the axial current, the decomposition is into longintudinal and transverse 
tensor structures. 

Subtracted dispersion relations can be derived for the scalar amplitude T^^(g^) assuming 
to be an analytic function of in the entire complex plane except for possible poles on and a cut 
along the positive, real axis starting at some threshold, see Fig. [T] applied to T^^ ijp') /[{p^)'^ {p^ — g^)], 
{N > 0), and often to a circular origin-concentric contour of infinite radius deformed to spare out the 
positive real axis. Closing the circle at a finite radius Sq < oo leads to so-called finite-energy sum rules 
|b2|, ES] which we will not focus on here. We have 



f_n2W ^oo Im T^^r si d^(T^^(t) xt)/dt 

where = —q^ > 0, > and it has been assumed that the integral along the circle vanishes, i.e. for 
a given the amplitude T^^ has an according behavior at infinity. For = the second term on the 

right-hand side of Eq. (jlj) is omitted. The quantity ^ is called subtraction constant. It 

is irrelevant if a so-called Borel transformation is applied to the dispersion relation Eq. (@]), see below. 

The discontinuity of T^^{q^) across the cut and the residues of isolated poles are related to the 
hadronic spectral function p^^{s) which follows from an insertion of a complete set of hadron states 
inbetween the currents of Eq. (^. It is proportional to the total hadronic cross section in the considered 
channel. This is also known under the name Optical Theorem. For example, the total hadronic cross 
section cr for the process 

e~^e~ ^ 7 — > hadrons with isospin 1=1 (5) 
according to the Optical Theorem is related to Im T^^ (associated with the current of Eq. (0)) as 

a{e'^e ^ 7 — hadrons with isospin 1=1) = ^Im T^^(s) , (6) 
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where a denotes the electromagnetic fine-structure constant. In sum rule applications one either sub- 
stitutes a fit to the low-energy spectrum together with the perturbative QCD spectrum for s greater 
than the threshold Sq for Im T^^{s) oc p^^{s) to determine unknown vacuum parameters, see below, or 
one assumes the dominance of the lowest-lying hadron within the resonance region in a given channel 
and determines its properties from known vacuum parameters. This leads us to the SVZ approach to 
the left-hand side of Eq. (jU employing the so-called Operator Product Expansion (OPE) in QCD. 

The OPE was originally proposed by K. Wilson |^ in the framework of a so-called skeleton theory 
for hadrons. The claim is that for a given theory a nonlocal product of composite operators Oi{x)O2{0) 
can usefully be expanded into a series 

Oi(x)O2(0)= Ec-(a:)O-(0). (7) 

D,iD 

involving local operators O]^ (0) of increasing mass- dimension D and c-number coefficients (x) - the 
so-called Wilson coefficients - as long as the (Euclidean) distance |x| is small compared to the inverse 
of the highest dynamical mass scale in the theory. The OPE can be shown to exist in perturbation 
theory, see for example [3H]. In Eq. ((7)) the index runs over all possible, independent operators of 
dimension D. Taking the vacuum average of Eq. ((7)), only scalar contributions survive as a consequence 
of the Poincare invariance of this state. 

Applying the OPE to a correlator of gauge invariant currents in QCD, the scalar amplitudes in a 
given current correlator are expanded into a series of the form Eq. ((7)). The participating operators, 
such as 1, f G^j^Gf^'' and m{uu + dd), are gauge invariant, and the associated Wilson coefficients are 
perturbatively expanded in powers of a^. By definition, the expectation in the perturbative vacuum of 
the expansion singles out the unit operator 1. 

In general, Wilson coefficients and operator averages depend on a normalization scale /i at which 
the perturbatively calculable short-distance effects contained in the former are separated from the 
parametrized long-distance behavior induced by nonperturbative fluctuations in the latter. If the cur- 
rents of interest are conserved then the associated correlator should not depend on any factorization 
convention. In this case a residual /z dependence of the OPE is an artefact of the truncation of the 
as series in each coefficient function. An improved scale dependence of the Wilson coefficients can 
be obtained by solving perturbative renormalization group equations. At higher order in opera- 
tors mix and anomalous dimension matrices must be diagonalized to decouple the system of evolution 
equations ^|, see also jHSl- Vacuum averages of nontrivial QCD operators defined at /z are universal, 
that is, channel-independent parameters. They are the QCD condensates. These parameters have to 
be extracted from experiment, and they induce (up to logarithm's arising from anomalous operator 
dimensions) (negative) power corrections in Q. Already in the seminal paper [2Tj it was realized that 
the OPE in QCD can at best be an asymptotic expansion since small-sized instantons spoil the naive 

power counting for operators of the form G^^^ starting at the critical dimension D = 12. 

This estimate was obtained from a dilute gas approximation for the instanton ensemble. Ever since, 
the discussion about the validity of the OPE in the light of its possible asymptotic nature, nonpertur- 
bative, nonlocal effects and the tightly related issue of local quark hadron duality has never faded away 
(see for example [Ml 1112 1113 Ell ESI EHl EH 112 iSl). No definite conclusions have been reached 
although valuable insights were obtained in model approaches. To shed light on this discussion is in 
the main purpose of the present review. It is stressed here, however, that the phenomenological success 
of many QCD sum-rule applications using 'practical' OPEs, typically truncated at D = 6, 8, supports 
the pragmatic approach of Shifman, Vainshtein and Zakharov of 'scratching' into the nonperturbative 
regime by ignoring questions of OPE convergence. For example, the prediction of the mass of the rjc 
at 3.0 GeV by a sum rule in the pseudoscalar cc channel matches almost precisely the result of a later 
experiment [¥T] . 
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Let us now discuss how Wilson coefficients are computed in practice. The perturbative part of a 
mesonic current correlator is the usual vacuum polarization function. In the case of electromagnetic 
currents this function was calculated to 4- loop accuracy in [1331 For the calculation of the Wilson 
coefficients appearing in power corrections two methods are known. The so-called background-field 
method |13Uj uses a Wick expansion of the time-ordered product of currents and interaction Lagrangians. 
Quark propagation is considered in a gauge-field background and the Fock-Schwinger or fixed-point 
gauge x^Af^{x) = is used. In the background-field method the gauge field and the quark field away from 
the origin are expressed in powers of adjoint and fundamental covariant derivatives acting on the field 
strength and the quark field at the origin, respectively. It is tacitly assumed in all practical applications 
that the gauge invariance of nonlocal operator products, as they arise in the Wick expansion, is ensured 
by straight Wilson lines 



which connect a field at Zi with a field at the base point zero (see |14Uj for an exhaustive discussion). 
The symbol V in Eq. (jH)) denotes the path-ordering prescription. The gauge invariance of the resulting 
nonlocal operators is at the price of allowing for gauge parallel transport only along straight lines 
connecting to the base point zero. The method does not consider fiuctuating gluons and therefore the 
computation of radiative corrections in Wilson coefficients is out of reach. It has been applied to the 
computation of Wilson coefficients for the gluonic operators of the form g^G^, g^G^, and g^{DGY = g^j"^ 
which induce the relevant power corrections in heavy quark correlators |137| I138t ll4Uj . A second 
method, which allows for the consideration of radiative corrections, works as follows. Project out a 
particular Wilson coefficient in an OPE of interest by sandwiching the associated T product of currents 
with (hypothetical) external quark and/or gluon states, separate off the structure belonging to the 
corresponding operator, and take the limit of vanishing external momenta. The calculation can be 
carried out at any order in perturbation theory - radiative corrections can be considered. This method 
was used in and in |1351 1136j together with the background method. 

Having discussed the two sides of the dispersion relation (@)) we will now review refinement procedures 
for the practical evaluation of sum rules. Typically, precise information about the spectral function in 
a given channel is limited to the lowest lying resonance. It can thus be important to apply a weight 
function to the spectral function which suppresses the high-energy tail to sufficiently reduce the error 
in the spectral integral. A Borel transformation, which acts on the sum rule ^ by an application of 
the operator 



is designed to do precisely this. On the spectral side it generates an exponential ~ exp[— s/M^] as 
opposed to the power suppression ~ l/(s -|- Q'^) at high energy and on the OPE side power corrections 
in the Borel parameter M are factorially suppressed in the mass-dimension. A physical parameter, which 
is extracted from the sum rule, should be independent of M within the so-called stability window. This 
region of values for M is interpreted as the interval where the spectral side of the sum rule indeed 
approximates the OPE expression. Hadron parameters are extracted in this region. 

It sometimes proves useful to suppress the high-energy tail of the spectral function by a higher 
integer power of s and then look at the insensitivity of the sum rule under variations in this power. In 
this case one considers various derivatives^ of the sum rule Q with respect to the external momentum 
Q^. This leads to a set of moment sum rules. 

^In contrast to T^^{Q^) the first derivative dQ2T^^ (Q^), the Adler function, needs no subtraction of the divergence 
related to the fermion loop at 0{a^). 








(9) 
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4 Seminal sum-rule examples 



Violations of local quark-liadron duality are believed to have a mild effect on QCD sum rules. We 
review specific QCD sum rules for (axial)vector mesons in this section. In view of future and ongoing 
experiments at facilities like the large hadron collider (LHC) and the relativistic heavy ion collider 
(RHIC), respectively, we also discuss the in-medium approach to QCD sum rules for (axial) vector 
current-current correlators. 

4.1 Borel sum rules for light (axial) vector meson currents 

We consider the following currents, which interpolate to the p, co", and Ai resonances: 

Jm'' = ^(^7/.w T dl^,d) , j^' = TV^{u-ff,j5U - d-i^,1->d) . (10) 

The operation TV projects the axial current onto its conserved part which allows to discard the pion 
pole contribution in the spectral function of its correlator. After separating off a transverse tensor 
structure q^qy — q^g^y the OPEs for the current- current correlators T^'^ , and T^^ in vacuum are given 



as 



(m7/.75A"m T c?7^75A"c/) ) (jj) - ((u-ff^X^u + d-f^X'^d) ^ qit^X'^q) (jj.) 



na. 



q=u,d,s 



2Q^ 

In^ y ' n I 2Q* ^'^^ ' ' 2AQ'^ \ n ^ ^^^^ 



q=u,d,s 



where a one- loop correction is included in the perturbative part, the Wilson coefficients in the power 
corrections are given without radiative corrections, fi denotes the normalization point, degenerate light- 
quark masses mq = rriu = md have been assumed, and A° are SU(3) generators in the fundamental 
representation normalized to trA°A* = 25°'^. In contrast to dimension four the operator averages at 
D = Q depend on the normalization point /x logarithmically |2II. For simplicity and because effective 
anomalous operator dimensions are small, we neglect this dependence here. 

On the spectral side, we use a zero-width model for the lowest resonance and approximate higher 
resonances by a (channel-dependent) perturbative continuum starting at threshold s™^^ (see for example 
|34j). The resonance part involves a meson-to- vacuum element. In the case of vector mesons it is 
parametrized as 

(0|j;7'"|mes(p)) = fmesmmese^, (12) 

where fmes {fries = p, u) is a decay constant and the meson's polarization vector. The contributions 
of the pion continuum to the spectral functions for s < sq are for practical reasons neglected here^. 
In Fig. 121 the experimental data for the sum of p and Ai spectral functions is shown for center-of- 
mass energies up to y/s = a/SGcV. The two resonance peaks at ~ (0.77)^ GeV^ = 0.59 GeV^ and 
~ (1.26 GeV)^ = 1.59 GeV^ are clearly visible. A similar situation holds in the u channel. After 
Borel transformation we obtain the following sum rules 

^At finite temperatures T this contributions becomes important for T > m^. 
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Figure 2: The sum of p and Ai spectral functions as measured by the Aleph collaboration 
in r decays. Plot taken from |179j . 



^ ^ M2 



(13) 



where the left-hand side of Eq. fjl3p can is represented by the Borel transforms of the respective OPEs 
(applying the operator ^ term by term). To extract information on the meson parameters one solves 
the sum rule, Eq. (fT3| . for the resonance piece R = fmes^~'^"^'"'^^^ ■ Assuming the coupling constants 
fmes and the resonance masses rrimes to be sufficiently insensitive^ to changes in the Borel parameter 
T = one can solve for rrimes by performing a logarithmic derivative 

d 

"^Ls = -^log^- (14) 
Using the following numerical values for the condensates 

{qq) {i2 = 1 GeV) = -(250 MeV)^ 

(yu = 1 GeV) = 0.012 GeV\ (15) 



TC 

a continuum threshold Sq = 1.5 GeV^, and approximating the four-quark condensate by squares of 
chiral condensate by assuming exact vacuum saturation (taking place in the limit of a large number 
of colors, Nc oo), one obtains a dependence of rup on r as shown in Fig. El We have neglected 
the contribution of the quark condensate at D=4 since it is numerically small compared to the gluon- 
condensate rriq/AQCD -C 1. Fig. 01 clearly indicates that the sum rule for the p meson is stable in the 
Borel mass for values around = 0.8 GeV^. One obtains qualitatively similar results in the u and Ai 
channels |22j . 

•^This is a a bootstrap-like approach to the sum rule. We first assume insensitivity of hadron parameters to variations 
in the Borel parameter and show subsequently that this assumption is self-consistent. For rUp this is shown in Fig.O The 
value of So below is chosen such that the size of the window of insensitivity is maximized. 

"^The quark condensate is determined by the pion decay constant /^r, the pion mass m.^-, the light current-quark mass 
by virtue of the Gell-Mann-Oakes-Renner relation and the gluon condensate can be extracted from ratios-of- moments 
in the J/\I' channel, see and next section. 
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Figure 3: The dependence of rup on the Borel parameter r. 



4.2 Quarkonia moment sum rules 

So far we have considered hght-quark channels using Borel transformed sum rules. In the case of mesons 
containing a heavy-quark pair moment sum rules turn out to be useful [HI]. The methods of Ref. [HI] 
have been developed further over the years and often applied to determine the mass of the bottom 
quark (for a review see [^). To be sensitive to the value of the quark mass the focus in the past was 
mainly on large moments since these suppress the perturbative continuum I139| I14()j . In general, 
moment sum rules require a precise analysis of the quarkonium threshold and an according definition 
of the quark mass (TO]. We briefly introduce the method here since we will rely on it in Sec. l7.3.!Zl 

One starts with the correlator where now light-quark currents of definite SU(3)i? quantum 
numbers are replaced by heavy-quark currents such as = h'^^h, h denoting one of the heavy-quark 
fields c, 6, t. Since this current is conserved we may again separate off a transverse structure {qfj.qu—q'^gfj.u) 
and then only consider the scalar amplitude T^^. In contrast to the light-quark case, where the OPE 



essentially is an expansion in 



^QCD 
Q 



the scale that power-suppresses nonperturbative corrections is 



naturally given by the heavy-quark mass m^. One therefore expands both sides of the sum rule in 
powers of the dimensionless parameter z = (5^/(4m|), 



(16) 



where Qc denotes the electric charge of the heavy quark, and compares coefficients C„. These can be 
expressed in terms of the moments J^n, defined as 



as follows 



-Mn = -r 



Mn 




(17) 



Q2=o 



Cn ■ 



The perturbative part of the coefficient C„ is known up to order and n = 8 \25b[ Wl\ IHHI • It can be 
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written as 



Cn = + ^ + Ci-) W) + + CT'l... + Ci-)4J , (19) 

where Im^ is a short-hand for \og{m\/ fi^) and the coefficients C^"^ are hsted up to = 8 in [69j. The 
nonperturbative part of the moment induced by the condensates (s'^G" ) , (g^ f"'^'^G'tyG\,)fl'^ ^ 



and (^g'^j'^j^'j, j/i being the hght-flavor singlet current, was calculated in |138t ll4Uj by using the 
background-field method. In our convention, Eq. (jl7p . the corresponding expressions are listed in |73j . 
Low-n moments are more sensitive to the spectral continuum, large-n moments to the resonance part 
of the spectrum. The former probe the relativistic part of the spectrum, and an expansion in as is 
appropriate. The latter probe the nonrelativistic physics in the quarkonium threshold region. The 
expansion parameter is modified by the velocity of the heavy quark and given as ^JnoLs- For, say, 
n > 4 a resummation of the spectrum to all orders in \fnas by a Schroedinger equation for the bound 
state should be performed (for a review see |f7Tj). 

On the phenomenological side the correlator is expressed in terms of a dispersion integral. In terms 
of the heavy-quark pair cross section cTe+e-^cc+x and muon pair cross section cre+e--*/^+/^- in e^e~ 
annihilations this leads to the following expression for the moments 

127r2Q2 J ,n+i a,^^__^^^^_(s) ' ^ ' 

It is obvious from Eq. ()2()|1 that with increasing n the hadron spectrum is probed at lower and lower 
center-of-mass energy y/s. To reliably extract the quark-mass parameter nih at high n therefore needs 
a high-precision treatment of the spectral function close to the threshold of quarkonium production. 
When rric was ffist estimated on parton level only the lowest four moments were considered |22j. 

If instead of analyzing the parameter rrih by using low moments one wants to estimate power correc- 
tions it is better to consider the ratio r„ = m^'\ adjacent moments since, contrary to the moments 
themselves, this quantity does not contain high powers of rrih for large n. Large n are needed to be 
sensitive to the nonperturbative information of the resonance region. The ratio-of-moment analysis 
thus is less sensitive to the error in the extraction of ruh and therefore more reliable. Historically, the 
value of the gluon condensate i^^G^^^G"^^ ~ 0.012 GeV^ was ffist estimated using a stability analysis 
of r„ for n = 2 ■ ■ ■ 10 in the ^/if) channel including resonances up to ?/'(4400) and the usual perturbative 
continuum model with threshold Sq = (4.2 GeV)^. In Sec. l7.3.51 we will have to go beyond the dimension 
four power correction when we analyze the 'running' of the gluon condensate in the framework of a 
reshuffied OPE containing nonlocal nonperturbative information. 



4.3 QCD sum rules in the (axial)vector meson channel at T, /i^ > 

In this section we address the case of a QCD current correlator in the (axial)vector-meson channel and its 
OPE inside a hot or dense medium. On the one hand, the high-temperature situation is of relevance since 
the properties of mesonic resonances (width and peak positions) in such an environment are measured 
in ongoing and future experiments at RHIC and LHC, respectively. To do this, the measurement of the, 
almost unperturbed by the nuclear environment, invariant mass spectrum of dileptons is carried out. 
Dileptons, produced in the early stages of a relativistic heavy-ion collision by vector meson decay, tell us 
about the mass and the width of the primary. A sudden 'melting' of the spectrum should unambiguously 
indentify the deconfinement transition. QCD sum rules are predestined to make predictions of the T 
dependence of spectral parameters |711[7g|7Hl|7gEIllHllEIlEHllH3Eniinil- On the other hand, 
there are interesting effects inside a cold baryon-rich nuclear medium which can be probed by heavy 
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ion colliders operating at lower center-of-mass energy. For example, isospin resonance mixing can be 
induced by an isospin-asymmetric environment. This is possibly accessible to a theoretical treatment 
using QCD sum rules and assuming the nuclear environment to be modelled by a dilute gas of nucleons 

The main focus of our brief review of the (axial)vector-meson sum rules in medium is on the oc- 
currence of a twist expansion^ when replacing the vacuum average over a bilocal current product by a 
Gibbs average 

(0| • ■ ■ |0) ((■ ■ ■)) = Z-iTre-^(^-^^«^) ■ ■ ■ , (21) 

where Z = ^xQ-P'^^-t^BQu) denotes the grand-canonical partition function for the QCD Hamiltonian 
H, fiB defines a baryon chemical potential, and Qb = J d^xj^ is the associated, conserved baryon 
charge. A resummation of such an expansion genuinely takes nonperturbative nonlocalities in the 
associated hadron states into account. At finite temperature T a small expansion parameter T/Q exists 
(T < Aqcd,Q ~ IGeV), and such a resummation is apparently not needed. The situation is quite 
different at finite nucleon density, where the expansion parameter is rrtp/Q {rrtp the proton mass) and a 
(partial) resummation of the twist expansion is imperative. The present section serves as a prerequiste 
to Seed where it is pointed out that the consideration of nonperturbative nonlocalities in vacuum 
OPEs is necessary for a good description of certain hadron properties. 

4.4 Light (axial)vector mesons at T > 

The pioneering work in formulating QCD sum rules for the case of finite temperature and/or baryon 
density was launched by A. I. Bochkarev and M. E. Shaposhnikov in 1985 [73]. In particular, the 
p-meson channel at finite temperature, as it was first treated in this paper, has been revisited several 
times over the years [7g|7g|7ZlEIll^lHnilZ3IH31IHIllHHllHainillBaE3 11021 because a number of 
substantial points were overlooked in [Tl] on both sides of the sum rule. 

Let us first discuss the basic points of the approach in [73] and its application to the p° channel. 
Instead of using a causal T product the formulation relies on a the retarded ordering of currents because 
of its more adequate analytical properties. In either case, the correlator of conserved currents T^i, can 
be decomposed into two invariants Ti(gQ = {u^q^Y ,q^) and T2(gg = {Ufj,q^)'^,q'^) since through the 
presence of a preferred rest frame a new covariant ^^=(1, 0,0,0) - the four-velocity of the heat bath - 
exists besides the external momentum q^, and therefore 

T^u{,q,u) = {q^qu - q^ g^,y)Ti + {u^ - uj'^){uu - ^^^)^2 , (22) 

where u = u^q^j,. In the limit g* — 0, however, Ti and T2 depend on each other, and it suffices to 
consider only one of them. By truncating the trace in the Gibbs average (PT|) to the vacuum and 
one-particle pion states (l-vr states) (dilute pion-gas approximation), for T < 160 MeV this is well 
justified due to the chiral gap in the spectrum (m,r ~ 140 MeV, nip ~ 770 MeV), a Borel sum rule for 
the channel at finite temperature was derived in [Z3] in the limit g — > 0. Evaluating the spectral 
side in this approximation, there appears a zero-temperature 2-7r continuum weighted by a function 
tanh(A/s/4T), usually neglected in vacuum sum rules, and, in addition, a scattering contribution to the 
spectral density, 

(x5{s)e{uj'^ -Aml)nB{uj/2T), (23) 

accounting for a l-vr intermediate state scattering off the current into a heat-bath pion and vice versa 
[76] . The p-meson contribution to the spectral function is as in vacuum since the l-7r-to-p matrix 
element of the current vanishes in the Gibbs-trace. At s > sq the spectrum is approximated by thermal 

^defined as mass dimension minus Lorentz spin 
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QCD perturbation theory. There is a quark-antiquark continuum and a scattering term. The latter 
arises from the interaction of the current with quarks in the heat bath. At T ~ 130 MeV it amounts to 
about three times the corresponding pion contribution. This shows the relative suppression of the pion 
scattering term in the hadronic part of the spectrum. 

On the OPE side, the T dependence of the dimension-six (four-quark) operators in Eq. (Illj) that 
contribute to the vacuum correlator was treated in [71] by applying the fluctuation-dissipation theorem, 
that is, by expressing the Gibbs averages of local operators through spectral functions in the respective 
channels 

((A(O)S(O))) = TT / ^coth(c./2T)p(c.,p,r), (24) 

where p is the absorptive part of the retarded correlator of the two currents A and B. In order to 
make sense of such a spectral function the original current product has to be Fierz rearranged into 
products of gauge invariant currents. Each of these products can then be evaluated using (j^ . For 
alternative ways of calculating pion averages over four-quark operators see below jTS] and [119J. For 
dimension-four operators this does not apply and some not too solid arguments were used to conclude 
that the T dependent part of their Gibbs averages is small [71] and therefore can be neglected. As a 
result of their Borel analysis Bochkarev and Shaposhnikov find that the both the p mass parameter and 
the spectral continuum threshold sq experience a drastic drop at T ~ 150 MeV. 

Three comments are in order: (i) On the spectral side in [71] the fact was not taken into account 
that at finite T the axial and the vector channel mix, and consequently that both resonances p and Ai 
should appear in the spectral function of the p° channel. In the dilute pion-gas approximation of [71] 
this was first shown in [TH] up to first order in the parameter e = Powers of e arise by reducing I-tt, 
2-7r, 3-7r, ■ ■ ■ states arising in the Gibbs-trace (j2T]l by means of the LSZ reduction formula, neglecting 
their momenta, and using PCAC and current algebra in the chiral limit. Taking into account only these 
finite-T corrections, the correlators T^^^^ (g, T) are expressible as a superposition of vector and axial- 
vector correlators T^i^^^q, 0) at zero temperature. Up to order only the e expansion contributes, and 
one obtains, 

T;,(g,T) = (l-6)T;,(g,0) + 6T;tn?,0), 

T^J{q,T) = (l-e)r;tHg,0)+eT;,(g,0). (25) 

Thus the resonance poles of the p and Ai mesons do not move as a function of T (and would not to 
any order in e if this was the only expansion for T corrections). The relative weight of the Ai meson 
in the spectral integral, however, increases with growing temperature. Interestingly, it was found in 
[7E] by a Borel sum-rule analysis applied to a single-resonance spectral function (erroneously, since the 
p mass does not shift at order e) that the p mass shifts towards higher values as T grows. One still 
can interprete this result in terms of the spectral weight moving towards higher invariant mass-squared 
s as T increases (growing importance of the Ai resonance, see Eq. (|25p). Except for [84j none of the 
QCD sum rule analysis subsequently performed, which all assumes a single resonance in the spectrum 
of the (axial)vector channel, has reproduced this behavior. Going to order T'^, there is an order 
correction, arising from the zero-momentum 2-7r states, but also a correction of order (T^/Q^) which 
originates from finite pion momenta in the l-vr state. The latter was estimated in [H2] by expressing 
the two invariants T^2y defined in analogy to Eq. with — > and u; — > i/ = p'^q^ in the l-vr 
matrix element of the current product, in terms of the measured pion structure functions Fi 2{x,q'^), 
X = Q'^/2u, using dispersion relations. There is a (T^/Q^)^ correction in both the Lorentz invariant 
and violating parts Ti 2, respectively. These terms are induced by nonscalar condensates in the OPE 
leading us to the next comment on |74j . 

(ii) In [71] only scalar operators were considered in the OPE. The 0(4)-invariance is, however, 
reduced to an 0(3)- or rotational invariance by the presence of the heat bath, and thus a number 
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of additional operators are allowed to contribute to the OPE. This was first noticed in [7^] where a 
systematic twist-expansion was used to identify the relevant, nonscalar operators. The gluonic stress, 
contributing at dimension four, was further investigated in [S2] and in [23 |HH1 in view of operator 
mixing under a change of the renormalization scale. In |79^ the l-vr matrix elements of operators, such 
as parts of 6^00 at dimension four (6'^,^ denotes the QCD energy-momentum tensor) and operators of 
the type q'^oD^q at D = 6, with non-zero twist - all other operator averages were omitted because of 
non-calculability - were evaluated using pionic parton distribution functions in the leading-order scheme 
and the parametrization of jHSI at the sum rule scale = 1 GeV. As noticed above on general grounds, 
the thermal phase-space integrals over l-vr matrix elements of pure-quark scalar operators in the OPE 
are expressible in terms of zero-temperature condensates and expanded in powers of e only by applying 
the LSZ reduction formula, PCAC, and current algebra [7^]. The l-vr matrix elements of the operator 
Os/vrG^ can be calculated by using the QCD trace-anomaly [031, 0^ = — l/87r(ll — 2/3Np)asG'^ + 
Y^qTnqqq. As a result, the matrix element is proportional to the pion mass and thus vanishes in the 
chiral limit. Even for realistic pion masses the T-induced shift of the gluon condensate is negligible |79j . 
Assuming a single, narrow resonance plus scattering term plus continuum model for the spectral side, 
the Borel analysis of [7H] indicates a drastic decrease of the p mass and the continuum threshold sq at 
T ~ 160 - 170 MeV. Notice that with the choice ^so(T = 0) ~ 1.3GeV> mA,{T = 0) = 1.26 GeV in 
[TH] the Ai resonance can effectively be viewed as a part of the perturbative continuum. A mixing of 
the Ai and p channels was noticed in [ZO] on the OPE side, in accord with the general result in |76j . 
Similar results were obtained for the u and Ai channels. 

(iii) No attempt was made in [71] to consider a T dependence of the vacuum state itself. This 
is suggestive since parameters like /^r, entering the average over the pion state in Eq. ((211), have a T 
dependence which, in the case of is calculable in thermal chiral perturbation theory [7H]. Chiral 
perturbation theory also predicts the T dependence of the quark condensate [HHl- For an investigation 
of the Gell-Mann-Oakes-Renner relation at finite temperature and the calculation of the T dependence 
of the pion mass relying on finite-energy sum rules see [HH^. Nothing, however, is known about the T 
dependence of parton distribution functions. Assuming that a T dependence of the vacuum state is the 
dominating dependence of the matrix elements in the Gibbs average (and thus that the T dependence of 
l-TT matrix elements can be neglected - the effect would anyway be of higher order in T) and assuming 
that such a dependence arises only implicitly through a T dependence of the continuum threshold Sq 
(related to a T dependence of the QCD scale Aqcd), a scaling of the vacuum averages of operators 
with powers of their mass dimension 2n, (so(T)/so(0))", was introduced in |84- A T dependence of 
the vacuum state is, indeed, suggested by the condensation of center vortices in the confining phase of 
QCD. These topological objects can be viewed as coherent thermal states themselves, hence the (mild) 
T dependence of the ground state. Condensed center vortices are apparently responsible for quark 
confinement and chiral symmetry breaking The scaling with powers in so(^)/'So(0) introduced in 
|Hlj is an effective, phenomenological way of considering this effect. Note that this scaling does not 
capture the small effects of omitted, higher hadronic resonances in the Gibbs average. As a result, 
a positive "mass shift" similar to the one in jEl and, as a byproduct, a moderate drop of the gluon 
condensate around T = 160 MeV, which at least quahtatively is in accord with lattice results was 
obtained, compare Figs. HI and |S1 Notice that this decrease of about 30% (for so{T = 0) = 1.5 GeV^) as 
compared to the value at T = is practically entirely due to the vacuum average and not due to l-vr 
matrix elements. 

(iv) The usefulness of thermal, i.e. on-shell quarks (the scattering term), in [T^ is quite questionable 
at low temperatures. 

We have seen that thermal, practical OPEs of current-current correlators allow for additional, 0(3)- 
invariant, operators to appear. These operators arise in the Gibbs average from matrix elements over l-vr 
states with nonvanishing spatial momentum. An expansion in powers of T^/Q^ arises in the chiral limit. 
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Figure 4: The ratio of the T dependent to the vacuum gluon condensate as obtained in 
[Mj . The dot-dashed, sohd, and long-dashed curves correspond to continuum thresholds 
so(0) = 1.2 GeV^, so(0) = 1.5 GeV^, and so(0) = 1.8 GeV^, respectively. The dotted curve 
shows the effect of perturbative renormalization-group evolution on the vacuum part of four- 
quark operators when allowing for an effective normalization scale \J sq{T) / sq{Qi)Q. Taken 
from IHH. 



Temperature induced corrections in Gibbs averages over scalar operators can be organized as expansions 
in two parameters, T'^/6f^ and T'^/Q'^ (in the chiral limit). The former arises from the (repeated) use of 
the LSZ reduction formula, PCAC and current algebra treating the pion as a noninteracting, elementary 
particle in the soft limit; the latter arises from the structure of finite- momentum pions. Since practical 
vacuum OPEs are, roughly speaking, expansions in Aq^^/Q"^ we conclude that the "convergence" of 
the expansion is not threatened by finite-temperature effects. 

4.5 Vector mesons in nuclear matter 

The treatment of current correlation involving (axial)vector mesons in a cold and dense environment 
using QCD sum rules is technically analogous to the case of finite temperature. Much work has been 
devoted to the calculation of the change of the mass and width of light vector mesons in a baryon-rich 
environment [^3 EH EH 11001 IMl 11021 1103 COl 1103 COS COZI (for summaries see [M COS EH ) since 
these should be measurable in terms of the invariant-mass spectra of dileptons emitted in the course of 
a heavy-ion collision, for an analysis within the Walecka model see |112j . at a facility like SIS18 (GSI) 
with the HADES detector. For an adaption of vector-axialvector mixing to the situation of pions in a 
nuclear environment see |llUj . A sum- rule analysis of the p-uj mixing induced by an isospin asymmetric 
nucleon density Pn = Pp + Pn was performed in |113j . [58j. respectively. This effect occurs in vacuum 
due to the breaking of the SU(2)p symmetry by the different electric charges and masses of up- and 
down-quarks jll4j . For a sum-rule analysis of the off-shell situation see [OH]- It was shown in |113j 
that by an appropriate and realistic choice of the isospin asymmetry a„p of the nucleonic environment, 
defined as a„p = y^^, the vacuum mixing can either be compensated or enhanced. Most of the above- 
mentioned works are technically rather involved. A simple and beautiful discussion of the p meson 
(positive) mass shift in nuclei is, however, given in |12Uj . New developments concerning the treatment 
of nucleon matrix elements, in particular the ones of four-quark operators, deserve a review article in 
their own right, for a recent publication relying on the perturbative chiral quark model see |118j . Here 
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Figure 5: The T dependence of the gluon condensate as obtained on the lattice. Open circles 
and squares denote the result obtained for QCD with light and heavy djTiamical quarks. The 
dashed line shows the ideal-gas situation (a formula derived from a scaling argument in |94j). 
and the solid line is a fit to the lattice result for the T dependence of the gluon condensate. 
Both lines refer to SU(3) Yang-Mills theory. Taken from 



we focus on interesting OPE aspects at finite density which hint on a fundamental manifestation of 
strong interactions at purely Euclidean external momenta in both vacuum and hadron properties: the 
occurrence of strong nonperturbative correlations characterized by mass scales considerably larger than 
the perturbative scale Aqcd- Let us now briefly review some technical aspects of finite density sum 
rules. 

On the spectral side, a linear-density or dilute-gas approximation for the Gibbs average in ()21|) . 
which consists of taking into account only the 1-nucleon state besides the vacuum, again leads to the 
occurrence of a scattering term in the spectral function which is due to the scattering of a bath-nucleon 
off the current into an intermediate-state nucleon and vice versa. The question whether a treatment of 
in-medium resonance physics relying on the linear-density approximation is reliable for nucleon densities 
larger than the saturation density is open. Moreover, the consideration of finite vector-meson width in 
a pure sum-rule treatment of the resonance seems to be problematic |lU3j . The sum rule apparently 
contains too few information to predict both the density dependence of the resonance mass and the 
width. On the other hand, consistency of a spectral function calculated in the framework of an effective 
chiral theory with the in- medium OPE of the correlator of the associated currents was obtained in |l(Jlj . 

On the OPE side, 0(3)-invariant operators contribute and can be organized in a twist expansion. 
Their nucleon averages are expressed in terms of integrals over nucleonic quark parton distributions, and 
a new expansion parameter, qI-, emerges. In practice, one omits twist-four and also mixed operators 
due to the very limited information about their nucleon averages. The nucleon average over qq and 
^G"^ are determined by the nucleon a term and by using the QCD trace anomaly, respectively, see 
|lU8j . The treatment of nucleon averages over scalar four-quark operators is not as straight-forward as 
in the pionic case where chiral symmetry fixes these matrix elements in terms of vacuum averages. One 
way of proceeding is a mean-field like approximation^ (MFA) adjusted to the linear-density treatment 
llUUj . For a treatment beyond the linear-density approximation methods have been worked out 
in |1081 llllj . The status of the MFA is quite obscure (for a recent discussion see |106| 1107] where 

^We use the same terminology as in P^ . 
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the strong sensitivity of the in-medium mass-shifts of p and uj mesons on the value of the in-medium 
four-quark condensates is stressed). 

The evaluation of the Borel sum rules in the p channel yields a decrease of the p mass [011 IHZl 1101 HOOj 
with increasing density. As for a the behavior of width and mass no definite conclusion is possible 
|1U4[ llU5j . The old results for the uo channel in |^ llUm 1^ . where in comparison to the p channel an 
enhancement of the screening term by a factor of 9 was overlooked |1UH I113j and a negative shift of the 
resonance mass was obtained, are in clear contradiction to more recent analysis |113| llU7j which points 
towards a positive mass shift. Note, however, that the calculation of the u mass shift in |l(Jlj . which 
is based on a chiral, effective theory, also indicates a negative sign. Consistency with the OPE in this 
case was reached by applying nuclear ground-state saturation in the same way as in the vacuum: 

{n{p)\{q^,^,X''qy\n) = - {n{p)\{q^,X^qy\n) = ^K{p) {n{p)\qq)\nf . (26) 

This approximation is different from the MFA. In Eq. a density dependence of the correction factor 
K is allowed for. 

Let us make some summarizing comments on practical OPEs at finite nucleon density, (i) As we 
have seen, a new expansion parameter, arises in the Gibbs averages over finite- twist, nonscalar 
operators. Recalling that rrip ~ 940 MeV and that the external momentum Q (or the Borel parameter 
M) should be not much larger than ~ 1 GeV to be sensitive to resonance physics and associated power 
corrections in the OPE, we must conclude that a naive expansion is hardly controlled. However, as it 
was shown in |115j . a summation of the twist-two correction to all orders in ^ appears to resolve this 

problem. Such a successful, partial summation of powers of gf- stresses the need to take nonperturbative 
nonlocalities in nucleon matrix elements into account. That this is not only true for the nucleon or, 
more generally, for any sufficiently stable hadron will be shown in detail in Sec. [T] where nonlocalities in 
vacuum matrix elements are imperative for a good description of certain hadronic properties. A very 
thorough discussion of the limitations of a local expansion of current correlators in the framework of 
nucleonic sum rules at isospin-symmetric finite baryonic density and of possible ways of improvement 
is performed in This discussion rests on the pioneering work |117j on QCD sum rules for the 

nucleon at finite baryonic density, (ii) The screening term can dominate the density dependent part of 
the sum rule (for example in the u channel or for the mixed p-u correlator with mean-field treatment of 
four-quark operator averages) |113j . We can take this as a general indication that most of the density 
dependence of the resonance parameters is induced by the hadronic model for the rest of the spectral 
function and not by QCD parameters. So the situation is reversed as compared to vacuum sum-rules, 
where information on the lowest resonance is obtained in terms of QCD parameters and not in terms 
of extra hadronic information, (iii) The status of the mean-field treatment of nucleonic matrix over 
four-quark operators llUUj is unclear. It was shown in |107j how a change by a factor of four in the 
contribution of four-quark operators can already change the sign of the mass-shift of the u resonance 
at nuclear saturation density. 



5 OPE and Renormalons 

In renormalized perturbation theory the divergent large-order behavior in correlators like the one in 
(P) can be related to power corrections of these objects [HHl IHH tHHl ISSl EZl 1^ EHI- In this section we 
very briefiy discuss the origin of this phenomenon and applications in QCD. We strongly draw upon 
the review by Beneke which contains the relevant references up to the year 1999. We will explicitly 
refer to only some of the subsequent developments in applications of renormalons. 
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In a power-in-as perturbative expansion up to order of, say, a two-point current correlator ^ 

T{a.) = E (27) 

n=0 

certain classes of diagrams, which we assume to dominate the expansion in ag in QCD, are associated 
with factorially-in-77, increasing coefficients r„ ~ Ka"'n^n\, (a,b,K constants), at large n . In this case 
the expansion would be asymptotic, that is, there exists a truncation A^* < oo which minimizes the 
truncation error. In gauge theories like QCD no proof is available for this asymptotic behavior. 

To have a sensible definition of a divergent series with factorially growing coefficients r„ it is useful 
to first look at the Borel transform of this series. For the series in Eq. ()27p it is defined as 

oo j.n 

B[Tm = Y.rn-- (28) 

n=0 

For a B[T]{t), which has no non-integrable divergences on the positive, real t axis, and which does not 
increase too strongly for t — >■ oo, one can define the Borel integral as 

/•oo 

T{as) = / dte-'^"^B\Tm. (29) 
Jo 

If T(as) exists then it defines the Borel sum of the original series T{as). If B[T]{t) has poles, which 
would then be a map of the diverging behavior of the series T{as), in the domain t > then one can still 
define a Borel integral for B[T]{t) by deforming the integration path in the complex t plane such that 
these singularities are circumvented. As a result, the Borel sum T{as) usually acquires an imaginary 
part. There is, however, no unique deformation prescription - poles can be circumvented by deforming 
to positive or negative imaginary values of t - which could be obtained from first principles in QCD 
perturbation theory. The difference between the two possible prescriptions embodies an ambiguity of 
the Borel integral which generically can be removed by adding exponentially small terms ~ 
to the power series T{as). One refers to the poles on the real t axis, which originate form factorially 
diverging coefficients in the perturbative expansion, as renormalon poles. 

Following the presentation in [SUj let us now look more specifically at how such singularities arise. 
We consider the Adler function, which is defined as 

DiQ')=An'Q'^^T{Q'), (30) 

because it is free of divergences related to the outer fermion loop. In Eq. (|Hn|) T{Q^) is defined as in 
Eq.dSD. 

More specifically we are only interested in contributions arising from chains of fermion bubbles 
as in Fig.ini At each order in these contributions are gauge invariant by themselves. The QCD 
renormalized fermion bubble leads to the following fermion-bubble-chain induced expression for the 
Adler function 

oo „oo 

where ^ = — /c^/Q^, k denoting the momentum flowing through the chain. The fermionic contribution 
to the (scheme independent) one-loop QCD (5 function is deflned as /5o/ = > 0, ^ Q denotes the 

normalization point, and the internal fermion-loop subtraction has been performed in the MS scheme. 
The function F{C,) is known exactly. It implies that for large n the integrand in Eq. (|3ip is dominated 
by ^ < 1 and ^ > 1. In the former case F(0 ~ f^^ and the latter F(0 ~ ^^^^ (log^^ + §)• This 
^In what follows the nonexistence of a constant term in Eq. (|27|) is inessential. 



PofOis log ^' 



(31) 
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leads to the following approximate (the low-n contributions are not well approximated) expansion in 
a, of the Adler function 



Dm = - E 



a 



n+l on 



0/ 



n=0 



(-2)-"n! + [n + — 



(32) 



The first (sign alternating since /5gj > 0) and second (sign non-alternating since /^qj > 0) terms in the 
square brackets in Eq. (jH^ are due to the ^ <^ 1 and ^ 1 contributions to the integral in Eq. (|^. 
respectively. The Borel transform of Eq. ()32j) reads 



B[D]{i 



2 /g^_ 



5/3 



1 I 4g^^-5/3 

2 - 97r /i2 



1 + 



+ 



5 1 



61 + t; 



(33) 



where v == —(3oft. The pole at f = 2, which is related to the behavior at small chain momenta, 
^ -C 1, is called first infrared (IR) renormalon whereas the single and the double pole at f = — 1, which 
originated from large chain momenta, ^ ^ 1, is called first ultraviolet (UV) renormalon. According 
to Eq. ()29|) only the latter makes a contribution to the Borel integral and generates a negative linear 
power correction in Q~^. This is in contradiction to what we expect from the OPE where the leading 
power in is —2 arising in the chiral limit from the gluon condensate {^G'^^^G'^^^ . What went wrong? 
The problem can be traced back to the fact that in considering only (gauge invariant) fermionic bubble 
chains and consequently only looking at the fermion contribution to the full QCD (3 function we actually 
computed renormalon poles which are close to mimicking the large-n behavior of an Abelian theory. 
Working in a covariant gauge, one could naively add the gluon and ghost bubble chains. The result, 
however, would be gauge dependent. A gauge invariant prescription to incorporate non- Abelian effects 
into our large-order investigations is to simply replace /?o/ by the full one- loop coefficient l3o of the 
QCD j3 function. This prescription includes also non-bubble-chain diagrams. Since the sign of /?o is 
opposite to the one of /So/ the IR renormalon pole moves to the negative (or positive) v (or t) axis and 
thus contributes to the Borel integral whereas the UV renormalon pole ceases to make a contribution. 
As a result, the lowest nonperturbative correction is a power and induced by an IR renormalon. 
Replacing /5o/ by /5/ in Eq. (jHH) and performing the sum first yields 



de 



-5/6\ 



(34) 



where is the running coupling at one-loop. Thus the effect of a gauge invariant sum of diagrams 
including the fermion bubble-chain is to replace the chain by a single gluon line which couples to the 
external fermions via the one-loop running coupling as{^e~^^^)- Although our prescription Pqj — > Pf 
seems to be ad hoc it can be justified diagrammatically that renormalon poles are located at integer 
values of v (or values of t that are multiples of J^. Using Eq. it is easy to see that the first IR 
renormalon contribution to -D(Q^) is ~ Aq^j-^Q'^ with an ambiguous but /i-independent numerical factor 
where the scale Aqcd is a typical hadron scale. All this matches nicely with the OPE approach where 
lower power corrections are forbidden by the absence of the corresponding gauge invariant operators, 
and the operator f G^^.G^'^ is renormalization group invariant at one loop. One may then say that the 
first IR renormalon is factored into the condensate and is associated with chain momenta k ~ Aqcd ^ 
while the Borel summable UV renormalons, corresponding to momenta /c ~ Q S> /i, do contribute to 
the Wilson coefficients in an unambiguous way. 

As we have seen, the renormalon approach offers some insight into the structure of power corrections 
even though the coefficients of the power corrections are ambiguous. At present it is not clear whether 
the OPE is asymptotic or not, and IR renormalons can not be used to predict the convergence properties 
of the OPE as an expansion in powers of itself. In fact, they only indicate the very limited set 
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Figure 6: A bubble-chain diagram which contributes to the perturbative expansion of a 
two-point current correlator. 



of power corrections which are related to large-order perturbation theory. However, one may think of 
more possibilities for the generation of power corrections, namely, power corrections which are entirely 
beyond the reach of perturbation theory or power corrections arising in Wilson coefficients from short 
distances. From a comparison of an analytical continuation to < with experimental spectral 
functions it is obvious that the so-called practical OPE violates local quark-hadron duality in the sense 
that we have defined it in Sec.|21 The construction and phenomenological test of reorderings of the 
OPE, which contain summations of powers to all orders and yet allow for a factorization of the 
large momenta regimes as in ordinary OPEs, is discussed in Sec. 17. 31 

In the remainder of this section we list two modern phenomenological applications of renormalons. 
For event shape variables and fragmentation cross sections in lepton-pair annihilation into hadrons, 
which are not described by an OPE, the identification of power corrections is not clear cut. Resorting 
to the so-called large /Sq approximation for the perturbative expansion, power corrections to the log- 
arithmic scaling violations in these quantities were treated using renormalon resummations jHOI- For 
current correlators associated with lepton-pair annihilation and r decay the relative strength of power 
corrections in their respective OPEs can be predicted from the corresponding residues of IR renormalon 
poles in a given scheme assuming that renormalons are the only source for these corrections. Comparing 
the large /3o approximation, in which this program is carried out, with known, low-order exact results 
only provides a partial justification for this approximation. This approach is an interesting model for 
power corrections and provides semi-quantitative insights, for a excellent discussion of this issue see (SO] 
and references therein. 

6 Violation of local quark-hadron 

In this second part a review of the experimentally measured violation of local quark-hadron duality in 
inclusive processes is given. Attempts to understand this violation using the model of current correlation 
with quarks propagating in an instanton background are reviewed. Finally, we discuss the issue in the 
framework of the 't Hooft model in the limit Nc —>■ oo, where the model is exactly solvable, and also 
for large, but finite N^. 

6.1 Experimental facts and lattice results 

Despite the practical successes of the use of the OPE in the framework of QCD sum rules, recall the 
moment analysis in the J/ijj channel fSec. l4.2|) . The theoretical status of this expansion in general 
field theories has never reached a satisfactory level, see |^ E3] for a discussion of scalar field theories 
with unstable vacua. In fact, it was even claimed in llHj as a result of an analysis of the 2D 0(N) 
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nonlinear a model that no cancellation between IR renormalons in the perturbative part of the OPE of 
the propagator with IR renormalons present in the condensate part takes place. This means that the 
definition of local condensates is ambiguous. 

In 4D QCD there is not yet an analytical way to decide on the role of perturbative contributions to 
the vacuum condensates. Direct calculations of current or field correlators were performed in (suitable 
limits of) various field-theory models and compared with the OPE ^"^^ it was found that the 
amount of perturbative contribution varies from model to model. Pragmatically assuming that the 
local condensates in QCD are dominated by nonperturbative effects, as it is done in any sum-rule 
application of the OPE, a clarification of the nature of this expansion in negative powers of the external, 
Euclidean momentum Q is still needed. This is, in particular, pressing in applications where analytical 
continuations of the OPE to the Minkowskian signature are needed as we will see below. 

Let us gather some experimental evidence that the inclusive spectra, which correspond to certain 
current- current correlators, do deviate substantially from the analytical continuation of their practical 
OPEs in the resonance region. We will consider three processes, (i) e^e~ annihilation into hadrons, 
(ii) axial-vector mediated, r decays into hadrons, and (iii) width difference in the decays of Bs and Bg 
mesons. 

(i) : In Fig.Elthe (electromagnetic) vector-current induced spectrum of e^e^ annihilation, 

_ atot{e^e~ hadrons) 
a[e+e ji+ji ) 

into hadrons is shown for y/s up to 14 GeV. It is obvious from the figure that a spectrum calculated 
from a continuation of the practical OPE, R = —^lm{Q'^ = —s — iO), violates local quark-hadron 
duality considerably within the resonance regions since the contribution of power corrections arising 
from operators with anomalous dimensions alter the perturbative result in Fig. [7| only in a smooth way 
at finite ^/s. 

(ii) : The isovector-axialvector-induced spectrum (labelled by A~) of r decay into non-strange 
hadrons can be expressed as follows |179j 



Ra~{s) 



Q\Vud\^SEw B(t 




NA-ds 



(36) 



where \Vud\ = 0.9752 ± 0.0007 is the CKM matrix element. Sew = 1.0194 ±0.004 accounts for radiative 
electroweak corrections, denotes the normalized invariant mass-squared distribution, and ~ 

1.77 GeV is the mass of the r lepton. In Fig.|Slthe spectrum Ra- of r decay as measured at LEP by 
the ALEPH collaboration |179j is shown. Again, there is no way for a naive OPE continuation into the 
Minkowskian domain to generate the behavior of the spectrum around the Ai resonance. 

(iii): Experimental information on Bg-Bs mixing is not yet available, but it will be investigated by 
CDF in the near future |149j . A theoretical prediction resting on the assumption of local quark-hadron 
duality in the OPE approach exists at next-to-leading order in as |154j . Following |155j we briefly 
give some theoretical background on why Bg-Bg mixing can be a testing ground for the violation of 
local quark-hadron duality. Since Bg may mix with its antiparticle Bg the two mass eigenstates Bh,l, 
which are linear combinations of Bg and Bg, have different masses. Am = — Mi ^ 0, and different 
inclusive decay widths, AP = Th — Pl 7^ 0. Am and AP can be related to the dispersive (absorptive) 
part of the Bg-Bg mixing amplitude, M12 (ri2), as follows 

Am = 2|Mi2|, AP = 2|ri2|cos0. (37) 
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Figure 7: The spectrum of e~^e annihilation into hadrons up to -y/s = 14 GeV. The cross- 
hatched band are the results of perturbation theory which agree nicely with the data in the 
continuum regions 2 GeV < < 3.1 GeV, 4.6 GeV < < 9.1 GeV, and 2^/^ > 11.2 GeV. 
The resonances in the region y/s < 2 GeV and in the vicinity of the cc and bb thresholds are 
out the reach of perturbation theory. Taken from |141j . 



Practically, we have AF = 2|Fi2| since the CP violating phase is very small in the Standard Model. 
By means of the optical theorem the Standard-Model expression for AF is 



AF 



(3J 



where T^e// denotes the effective Hamiltonian mediating transition between Bg and Bg in the Standard 
model after the heavy vector bosons and Zq have been integrated out using renormalization-group 
improved perturbation theory. The operators appearing in the decomposition of 'He/f are normalized at 
a scale fii = 0{mh). The variant of the OPE, which is used to estimate the right-hand side of Eq. (|38|) . is 
an expansion in inverse powers of the 6-quark mass mi, - the so-called heavy quark expansion (HQE) 15^]. 
Notice that this expansion relies on the validity of the HQE along the discontinuity in the Minkowskian 
domain. An experimental detection of a sizable CP asymmetry in the decays of Bg and Bg would signal 
new physics. To test the Standard Model it is thus extremely important to have a good understanding 
of nonperturbative QCD effects and in particular of the validity of local quark-hadron duality in the 
use of HQE. To lowest order in 1/mb one has 
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Figure 8: Spectrum of axialvector-induced, r decay into hadrons. Taken from |il79j . 



In Eq. fl39|) Gp denotes the Fermi constant, F and Fs are the imaginary parts of the Wilson coefficients 
in the leading-order in as expansion, and Q{Qs) are the /S.B = 2 operators Sj7^(l ±75)6.jSj7'^(l ±75)6^-. 
The matrix elements of Q and Qs are parametrized as 

Bs\Q{f,,)\Bs) = -flMlm,) , lBs\Qs{^i2)\Bs) = --flMl^-——^^—-^Bs{^i2) (40) 



where /^^ and Mb^ are the decay constant and the mass of the Bg meson, and the masses mb{^2) 
and ms{ii2) are defined in the MS scheme, and the normalization scale 112 = 0{mb). In the vacuum- 
saturation approximation the "bag" factors B{fi2), Bs{fi2) are equal to one. 

At next-to-leading order in there are already seven operators^ which describe nonlocal contri- 
butions to the transition. Including corrections |157j . the final answer for the quantity AF/F, 
F = 1/2(Fl - Fh), reads [THH] 



F I 245 MeV / 



[(0.234 ± 0m5)Bsimb) - 0.080 ± 0.020] (41) 



where rribimi,) + msimb) = 4.3 GeV (MS scheme) and ml/ml = 0.085 have been used. Due to its 
tiny numerical value the contribution oc B has been neglected in Eq. Mlj) . With the result Jb^ = 
(245 ± 30) MeV (for a QCD sum rule determination see for example |152p of an unquenched lattice 
calculation |151j (two dynamical fermion flavors) and the result Bs{mb) = 0.87 ± 0.09 of a quenched 
lattice calculation |153j one obtains (lattice errors added linearly) 

AF , , 

— = 0.12 ±0.06. (42) 

In the limit of N^. ^ 00, where vacuum saturation is exact, and for Aqcd ^fe ~ <C nib one can 
show that local duality holds exactly |158j . In this case the result is ^ = 0.18 which is just in the 

^Dimensional regularization with anticommuting 75 matrices and the MS scheme is used (scheme dependence of the 
Wilson coefficients cancels against scheme dependence of the associated operator averages) in |154| . 
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upper error limit of Eq. ()42|) . It is thus clear that a future experimental detection {Nc = 3) of violations 
of local quark-hadron duality in the HQE for ^ and/or of New Physics needs a much more precise 
(lattice)determination of fs^B and f^_^Bs, see for an unquenched, Nf = 2 calculation of Bs |15Uj where 
the error (and the central value) are reduced in comparison to the result of the quenched calculation in 



6.2 Quark propagation in an instanton background 

Since no full, analytical solution of QCD exists, which would allow for a direct comparison of the OPE 
with the exact result (and make the OPE superfluous), one has to resort to models of the current-current 
correlator. It was proposed in [142j that in analogy to the cancellation of renormalon ambiguities, arising 
from a factorial growth of coefficients in the a<j expansion, by exponentially small terms e^^"'^ , the OPE 
may (at best) be asymptotic in the according expansion in ~ ^'^q° ■ To make sense of it, one would 
then have to add exponentially small terms of the form q^/^qcd ^hich possibly would cure the violation 
of local quark-hadron duality by practical OPEs. 

A seemingly reasonable approach to see whether there is some truth in this proposal is to consider 
the quark propagation |122j , inherent in a given correlator of currents with massless quarks, in a dilute- 
gas instanton-antiinstanton background |123| I124j or a general, (anti)selfdual background |127j (like a 
dilute gas of multiinstantons and multiantiinstantons). 

Let us briefly review the calculation of the correlator of electric currents in a dilute-gas instanton- 
antiinstanton background as it was performed in |124j putting the (anti)instanton into the singular 
gauge. In the dilute-gas approximation it is only necessary to regard a single quark flavor q of electric 
charge Q - a, sum over flavors can be performed at the end of the calculation. We consider the two-point 
correlator of the conserved current = Qq'j^q in Euclidean position space 

T,,{x) = {Tj^{x)jM) ■ (43) 

Considering, in a first step, quark-propagation in the background of a single (anti)instanton and disre- 
garding radiative corrections, the Lorentz-trace T{x) = Tl^{x) is simply given as 

T±(x,fi±) = -J2TTQlYSi{x,0,Qih^S^{0,x,Qi) (44) 

F 

where denotes the quark-propagator in the (anti)instanton background, the trace is over Dirac and 
color indices, the sum is over light quark fiavors, and Q± denotes the collective parameters of the 
(anti)instanton. The propagators are expanded in mass rrii around = 

a) = -MEliofe) + ^ liWMM + ,„ ^ liMllM ^ o(,„^) (45) 

where A is a nonvanishing eigenvalue of the Dirac operator i'^^D^, and the subscript '0' refers to the 
zero-mode contribution. It is important to keep the term linear in m when calculating T{x) in the limit 
m 0. The zero- mode part in Eq. is given in terms of 

= (|;)"' ((,_jr+,|)3/. {r,.i^'4(l±;.))^^^ (46) 

where x denotes a 4D unit vector, and a{t) is a Dirac index (fundamental SU(2) color index). The 
0(m°)-part in Eq. (jl^ . Sl{x,y,Q±), can be written as 



5° (x, y, n^) = YD',^±{x, y, ± 75) + (A±(x, y, fi±)D^7^)^(l t 75) (47) 
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where A±{x,y,fl±) denotes the propagator of a scalar, color-triplet particle, and {A±{x,y,Q±)Df^'-f'^) 
means that the covariant derivative acts from the right onto A±{x, y, This propagator is explicitly 
known \l'2'2\ . in singular gauge it reads 

/ plal{x-x±Yat{y-x±y\ ( pI Y^'^ , , 

1 + ^2, ^2 1 + i A2 (48) 

\ [x-x±)\y-x±Y )\ {y-x±yj 

where = {Rabcr'',^i)^, Rab G S0(3) is a (constant) rotation matrix in adjoint SU(2) color space, 
a'' ,{b = 1, 2, 3), denoting the Pauli matrices, p± the (anti)instanton radius, and x± is the center of the 
(anti)instanton. The 0(m) contribution in Eq. (jlHI), S]_{x,y,Q), can simply be expressed as 

Sl{x,y,Q)=m J d^z Sl{x, z,Q±)Si{z,y,Q±) . (49) 

Inserting the zero- mode expression ()46|) . and the zeroth- (first)- order in m expressions Eq. ()47p (Eq. ()49p) 
into Eq. ()44|) and only considering the part, which survives the limit m — 0, averaging over the color 
orientations of the instanton embedding into SU(3), subtracting the free current correlator Tq, perform- 
ing the integration over (anti) instanton centers and radii over the remainder, and taking into account 
the contribution from instantons and antiinstantons in this part, one arrives at the following expression 

[1211 

where -D(p) denotes the instanton density at one-loop perturbation theory (only gluonic fluctuations) 




(51) 



which can be interpreted as the number of instantons of size between p and p + dp per unit space-time 
volume. 

After separating off a factor -^q"^ (arising from the transverse tensor structure) in the Fourier 
transform of the entire correlator Tq + 6T{x,p) and after accounting for the Gaussian integration over 
fermionic fluctuations around the (anti)instanton we have 



T.Qf) Ug% + 16r:^ f^D{p)A{p) 

F \ P P 



1 1 y f 2Qp 

duK2 



(52) 



where A{p) denotes the associated (dimensionless) fermion determinant |125j . and K2 is a McDonald 
function. The integral over p is cut off at small p by the fermion determinant. For large p it is ill- 
behaved which signals that the dilute-gas approximation as well as the one-loop perturbative treatment 
of fluctuations breaks down. One usually introduces an upper cutoff pc ~ Y^cd by hand. 

Clearly, there is a dimension-four power correction arising from the first part of the integrand in 
Eq. (I52j) which can be associated with the gluon condensate. The second part of the integrand, however, 
is not a power correction. For asymptotic momenta Q — 00 it falls off exponentially and can be taken 
as an indication for the searched-for exponentially small terms needed to cure the OPE. The reader may 
wonder why there is only a dimension- four power correction in Eq. ()52|) . To answer this, let us recall 
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Figure 9: Experimental results for R{E) (right panel) and the result of Eq. obtained in 
the instanton model (left panel). The dashed line corresponds to the parton-model result. 
Plots taken from |142j . 



that the appearance of dimension-six condensates is associated with radiative corrections - the prefactor 
as before the four-quark operators refers to a gluon exchange initiated by the current-induced quarks. 
On our above treatment of quark propagation, however, we did only consider radiative corrections to 
the background-field but not to the quark propagation in this background. 

One may now continue Eq. ()52|) to the Minkowskian domain and determine its imaginary part to 
give a prediction for the ratio R = — -Im((5^ = — s — iO) in Eq. (j35|) . 

The result is 



R{E) = ^— (^R,(E) + R\E) 



1671= fdpDifi) 



' ' ' 2{Ep)^Jo Vv^ 
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(53) 



where E = ^/s, Ro{E) and R\E) denote the free particle and the instanton induced parts, respectively. 

Dip. 



and J2 is a Bessel function. In |142j the product -^r^A(p) was approximated by the simplest possible 



form 

^A(p)=rfoPo5(p-Po) (54) 

where a value po = 1.15 GeV~^ and do = 9 x 10~^ was adopted ^. A comparison of the function 
R{E) in ()54j) and the experimental results for R{E) is presented in FiglHl It is obvious that the part 
in Eq. ()53|) not contained in the practical OPE is responsible for the resonance-like behavior at low E. 
Although quantitatively the two plots in Fig.|nidiffer^° - after all it is clear that an incomplete dilute-gas 
approximation, recall the bold choice of the instanton weight in Eq. ()54j). is not a good approach - there 
is at least some qualitative agreement. 



6.3 Duality analysis in the 't Hooft model 

QCD in two dimensions (QCD2) considered in the limit Nc —>■ 00 with g'^Nc fixed - the so-called 't Hooft 
model |159j - is exactly solvable. At finite but large Nc a well controlled expansion in powers of 1/Nc 

^These numbers are obtained by requiring that the instanton induced contribution to the semileptonic width of D- 
meson decay are 50% of the parton-model prediction |142| . 

^"^Unfortunately, we have energy on the x-axis in the left panel and energy squared on the x-axis in the right panel. 
Even though this makes a direct comparison more cumbersome the author of the present review chose not to adapt the 
figures in ^142^ . 
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Figure 10: Planar diagrams contributing to the quark self-energy and the quark- ant iquark 
scattering amplitude. Thick lines denote gluon propagators (not dressed in the limit — > 
oo, g'^Nf. fixed) and thin, arrowed lines propagators refer to quark propagators. 



is available. For this reason it is the ideal testing ground for questions on local quark-hadron duality, 
namely, at large external momenta the practical OPE of some polarization operator (current correlator) 
can directly be compared with the asymptotically exact result, and duality violating contributions can 
be identified. A vast literature exists on the subject, see for example |16H I162| I163| ITMl 11651 H^^t H^^l 
I168t 117m I17H I172| I173t I174j . and not all contributions can be exphcitly referred to here. A review, 
which also discusses the string interpretation of the 't Hooft-model results, exists |175j . During the ten 
years or so the interest in the 't Hooft model was boosted by questions of duality- violations in the weak 
decay of heavy quark fiavors, see for example |167t ll7Ut 11691 ^TTTl 11701 ^^7^ ■ by the necessity to check the 
reliability of lattice calculations, see |1 73111 74j . by the need to estimate higher-twist corrections to parton 
distribution functions |165j . A dynamical understanding of chiral symmetry breaking in two dimensions 
was obtained relatively early |164j . In this section we are mainly concerned with duality violations in 
spectral functions based on the OPEs of current correlators when allowing for l/N^. corrections. 



6.3.1 Prerequisites 

Before going into the technical details of duality violations in the 't Hooft model we will here give a 
brief introduction into this model [159j. 

One considers the usual QCD Lagrangian 

C = -\gI,G^: + f{^l'D, - m^)q^ . (55) 

^ F 



Spacetime is two dimensional, the gauge group is \]{Nc) instead of SU(Ai'c), and the gauge coupling 
g has the dimension of a mass. Conveniently, one works in light-cone coordinates = i 2;°) 

and A± = ^ {Ai ± Aq) where = Xi and x° = — Xq. Imposing the (ghost-free) light-cone gauge 
A_ = y4+ = 0, one has = —d^A+ and Eq. reduces to 

C = \{d^Al)' + Y: fi^Yd, -mF + g^^A+)q^ . (56) 

^ F 

Taking as the new time direction, the field A°^ has no time-derivatives, and thus it is not dynamical. 
It will provide for a static Coulomb force between the quarks. Since 7^ = 7+ = 0, 7+7_ + 7-7+ = 2, 
and since the vertex in Eq. (j56j) comes with a 7_ one can eliminate the gamma matrices from the 
Feynman rules. Suppressing color indices, the gluon propagator is 1/fc^, the quark propagator is 
k-/{2k^k^ — m?p + ie), and the vertex is 2g. It was shown in |160j that in the limit ^ 00 with g'^Nf. 
fixed only planar diagrams with no fermion loops like the ones in Fig.llOlsurvive in the calculation of any 
amplitude. Due to this extreme simplification the equation for the quark self-energy (the rectangular 
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Figure 11: The Dyson series for the quark self-energy (upper figure) and the equation which 
determines it in the hmit Nc ^ oo, g'^Nc fixed. 




-X 



-X 



p-r 



Figure 12: Homogeneous Bethe-Salpeter equation which determines the spectrum of meson 
states in the hmit Nc — > oo, g'^N^ fixed. 



blob in Fig.llip can be written in untruncated form. 

To solve this equation requires, in intermediate steps, the introduction of a symmetric ultra-violet 
cutoff as well as an infra-red cutoff. The former is a consequence of the strong gauge fixing and has 
no physical interpretation. The removal of the latter in the final result, which does not depend on the 
ultra-violet cutoff, shifts the pole of the quark-propagator to infinity - one concludes that the spectrum 
has no single quark state. To determine the spectrum of quark-antiquark bound states one has to look 
at the homogeneous Bethe-Salpeter equation as depicted in Fig.lT^ 

Exploiting that the Coulomb force is instantaneous to separate the loop integrals and introducing 
the following dimensionless quantities (compare with Fig. lT^ 



where /i = [g'^N^/Tt and m„ is the mass of the nth meson, one obtains the 't Hooft equation for the 
mesonic wave functions 0„(x) 



with P denoting principle- value integration, P [l/(x — yY] = lime_»o \ — y + ie)^ + l/{x — y — «e)^]. 
In writing Eq. (j58j) we have assumed the masses of the participating quark and antiquark to be equal, 
mp = mp = m. It was shown in |159j that the "Hamiltonian" defined by the right-hand side of Eq. (j58|) 
is hermitian and positive definite (finite quark mass, 7 > 0) on the Hilbert space of functions which 
vanish at x=0,l like x^, (1 — x)^ where /3 is a root of vr/? cot('7r/?) = 1 — 7. The spectrum is discrete 
and is for large n approximated by 0„ ~ v^sin(7rnx) and /i^ = vr^n. For 7 = the lowest meson mass 
vanishes, the associated wave function is 0o = 1. 

Since the spectrum is real and positive definite resonances do not decay into one another - their 
width is zero. At large s (or n > no) the spectral function of a correlator T of currents, which couple 



2 2/2 2 2/2 



X = p-/r_ , 



(57) 




(58) 
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to all meson states equally, is therefore well approximated by equidistant, zero- width spikes: 



p{s) = ImT = const — - 6 — n] . (59) 

2 „>no Wf^ J 

6.3.2 Current-current correlator and spectral function at 0{1/Nc) 

Our discussion of local duality violation of current- current correlators in QCD2 beyond the limit Nc ^ 00 
relies on work |168j which uses the older results in |16H 11621 1163j . 

At finite Nc the quark- ant iquark bound states of the 't Hooft model are unstable. For a decay 
a 6 + c the width Ta at 0{1/Nc) is given as [Ml 11121 CESI 

- - ^ EE / ^""'^ (60) 



where I{ma, rrib, rric) = l/4[m^ — {mi) + mc)^][m^ — {rrih — mcY] and the meson coupling is given as 

g^bc = \l - (-1)-+-^+-] (/+^ + /-J . (61) 



N, 



c 



The parity of the ath meson is aa- The quantities are constants for on-shell decay. They are given 
by overlap integrals between the meson wave functions (f)a,b and the Bethe-Salpeter kernel, see |168j . 
For lack of better analytical knowledge Eq. was used in |168j with the asymptotic spectrum, even 
for a,b,c < uq. For small a one should not trust this approximation. For an investigation of duality 
violating components at large Q it is, however, justified. A numerical evaluation of Eq. ()60|) and a 
subsequent fit to a square-root dependence yields the following estimate 

r. = ^^^^^^^v/^ [1 + 0(1/.)]. (62) 

This knowledge of the dependence of width on n can be exploited to estimate the correlator T{q^) = 
i J d?xe^'^^ {Tj{x)j{0)) of the scalar current j = qq. In the limit iVc 00 it is given as 

00 2 

T{q') = - E 2 2^- • (63) 

Requiring this to match the leading perturbative order for — = — >■ cxd (duality in mesonic and 
quark description of asymptotic freedom in T(Q^ — > cxd)), one obtains |161j 

gl = N,T:p^, (nodd). (64) 

Inserting Eq. ()64|) into Eq. ()631) and taking the imaginary part, one arrives at an expression like in 
Eq. (j59|) . By means of a dispersion relation one can define ImT(g^ > 0) to T(g^) everywhere in the 
complex g^-plane, up to a constant. This gives |168j 



T(g2)-T(0) = --i^^(a), ^ = + - (65) 

277 27r^/i^ 2 



where V'(c") denotes the logarithmic derivative of Euler's Gamma function. So far we have discussed 
the narrow- width case [Nc —>■ 00). The finite-width case is treated in the Breit-Wigner approach by 
replacing — + ie ^ q"^ — m\ + S(g^) in Eq. Since ImE(g^ = m^) = m„r„ = one may 
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Figure 13: The spectral density ImT(s) in the 't Hooft model with the leading 1/iVc contri- 
bution to the resonance widths included. The center-of-mass momentum s is given in units 
of 27r^/x^. Plot taken from |,168] . 



take E(Q2) 
reads 



log 2^-^- As a result, the correlator with l/N^ corrections in the widths included 



r(g^)-T(0) 
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cr 



-A/iTr'^Nc)+l 



27rV 



+ 



(66) 



1 - A/{7T^N, 

Fig.lT!^ shows the spectral density ImT(s). The broadening of the resonances with increasing s is clearly 
visible. How does this result compare with the practical OPE? The OPE corresponding to Eq. (j66j] is 
known to be the following (asymptotic) expansion 



Nr. 



1 - A/{TT^Nr) 2n 



log a 



a 



n=l 



-ir-^- 



2n 



(67) 



The Bernoulli numbers behave at large n as Bn ~ {2n)\ which shows that the expansion is asymptotic. 



The variable cr can be obtained from a by replacing 



z. 



Consequently, the powers in Q in 



Eq. (jU7j) are slightly displaced at finite as compared to their integer values at Nc 
in this deviation a = _/lr , this introduces logarithmic corrections: 



2n~a 



l/Q^Y " {1 + alogQ'' + 



2n 



oo. Expanding 



(6J 



The logarithms in Eq. ()68|) lead to smooth contributions to the spectral density at Gnite s which, 
however, vanish in the limit Nc ^ oo. A direct (not using the asymptotic expansion ()67|) ) calculation 
of the spectral density reveals an oscillating component 2 exp[— as//x^] cos[s/(7r/i^)] for as S> /x^ which, 
however, is not suppressed by l/iV^. This contribution signals the violation of local quark-hadron 
duality. 



6.3.3 Decays of heavy mesons 

As a final testing ground for the study of the violation of local quark-hadron duality in the 't Hooft 
model we will now discuss the weak decay of heavy-light mesons. There has been a rather high research 
activity within the last few years investigating the issue of both global, see for example |167j , and local 
duality, see for example |17m I17H I169t I172j . This is natural since important lessons for the (3+1) 
dimensional case, which is targeted by present experiments at the B factories, can be drawn. 
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Figure 14: Lowest-order diagram for the transition amplitude in Eq. ()7()j) . The wave hne 
represents the propagation of the pseudoscalar massless particle which is composed of ipa 
and ip. Taken from |17Uj . 



Our discussion mainly follows |17()j . To describe weak decays the 't Hooft-model Lagrangian has 
to be supplemented by a weak-interaction part. Effectively, we can take it to be of the current-current 
form 

n^ff = ^ {qi,Q) {^ai'^,) (69) 

for two vector currents since the axial current reduces to the vector current in (1+1) dimensions. The 
constant G is the 2D analogue of the Fermi coupling in 4D. By means of the optical theorem the 
inclusive decay with reads 

= Im^ j d'x {HQ\TnXff{xW,ffmHQ) (70) 

where Hq denotes the state with the heavy meson at rest. Depending on whether we describe semilep- 
tonic or hadronic decay the fields ipafi in Eq. are either leptonic fields or quark fields. We only 
discuss the case where = 0. As was shown in |17Uj . the current ipal'^ipb can be substituted by 
e^'^ dv(l) / \/tt where denotes a pseudoscalar, massless field The corresponding lowest-order diagram 
for the transition amplitude is shown in Fig. El The parton-model result is 

Tq = ^ 71 

and can be obtained by replacing the heavy meson by the heavy quark to be able to calculate the average 
over the operator T° = CqqQQ with 2Imc^Q = Fq. Going beyond the parton-model approximation, 
the (spacetime integral of the) operator QQ needs to be expanded in powers of I/tuq. Up to ©(l/mg) 
one obtains |17()j 



Hq\QQ\Hq) ^ 1 {Hq\Q{-DI)Q\Hq) ^ {HQlQl.t-E.qrf^qlHQ) ^ ^^^^ 



2Mhq 2m^Q 2Mhq 2m^ 2Mhq \mf 



Including first-order in g"^ radiative corrections to the coupling Qq(j) simply amounts to a shift of the 
quark masses: mq^^ ^Q,q ~ /^^ where = limAr^^oo f^^c < co denotes the 't Hooft coupling. This 

9 9 

is actually true to all orders |17()j . As a consequence, we have to all orders that 2ImcQQ = ^ — | f- 
Under radiative corrections the expansion of cqqQQ up to ©(l/mg) thus is entirely due to the one 

^^In the leptonic as well as in the case where i/'ab £^re quark fields their polarization operator is exact at one-loop level 
owing to the facts that 7"7^7q — and that a product over an odd number of gamma matrices reduces to one. There is 
a pole at s = corresponding to a massless particle. In the case of quark fields it is known that in the 't Hooft model 
the only nonvanishing meson-current coupling of the vector current is the one to the massless meson. 
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Figure 15: The oscillatory part of Thq represented as the ratio Ar°*'^/(G^/3) versus the ratio 
Mfj /{n'^fS'^). Taken from \170l 



of the operator QQ, see Eq. ()72j) . As far as the first subleading, dimension two, four-fermion operator 
O^q = QTiQqT2q - Fi 2 denoting color and spinor matrices - is concerned^^ it can be shown that its 
Wilson coefficient c^g does not develop an imaginary part at leading order in g and at one loop. Thus 
this operator does not contribute to the total width up to 0(l/m^). There is a contribution of 0(l/m^) 
at two loop, and hence we have up to 0(l/m^) 



Eq. (|7Sj) can be compared with the following result of a calculation obtained by using the 't Hooft 
equation (now with mq and the mass of the spectator quark nisp 7^ ttiq) (170j 



Writing the operator QQ in terms of quark-components in the light-cone formalism and absorbing 
renormalization factors, the matrix element in Eq. ()73p can be expressed as a functional of (pHg such 
that Eq. ()74j) is reproduced. Up to ©(l/mg) the OPE prediction matches the expansion of the exact 
result. Also, the absence of a l/rrib correction, found by analyzing the OPE, can directly be verfied 
using the HQE of the 't Hooft equation in the approach of Ref. 

To derive Eq. ()74|) a sum over exclusive widths F„ for the decay into the nth meson had to be 
performed from which the off-shell part (m„ > Mhq) had to be subtracted afterwards. It can be 
estimated that the analytic part of the latter term is ©(l/mg). If there is a nonanalytic violation of 
local duality AF"'^'^ it must reside in the sum with m„ > Mhq ■ An estimate of AF"*'^ was made in |17(Jj 
for iTLsp < 13. In Fig.^1 AF°'^'^ is shown as a function of the heavy meson mass. The relative amplitude 
of the oscillatory part can be estimated as 



^^In contrast to (3+1) dimensional QCD, where the operator g^Qcr'^'-'^G^^Q is the first subleading operator in 2D it has 
dimension four and thus is not the first subleading operator. 




(73) 




(74) 




(75) 
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and hence it is strongly power-suppressed in the weak-couphng hmit Mhq ^ [3. 

To conclude, there is definitely an oscillatory, duality violating component in the inclusive decay 
width for the semileptonic or hadronic decay of a heavy meson in QCD2 considered in the limit iVc — 00, 
g'^Nc fixed. It is strongly power-suppressed for Mhq ^ P- 



7 OPE and nonperturbative nonlocality 

In the last section we used field-theory models, which are not too far from first-principle, realistic QCD, 
to pin down the terms absent in practical OPEs that lead to a violation of local quark-hadron dual- 
ity. This section approaches the problem from a more phenomenological side. By relating moments of 
hadronic light-cone wave functions (also called distribution amplitudes (DAs)) to the OPE represen- 
tation of an appropriate vacuum correlator, thereby gainig experimentally testable information on the 
wave functions in terms of perturbative QCD and condensates |184t 1185] . a phenomenological analysis 
of the validity of the OPE can be performed. Allowing for deviations from the locality of condensates - a 
possibility discussed since the early days of QCD sum riiles |TMI^ II!n MimiTM[Ti71[TH^ I7nil^ 
- strongly suggests that a truncated, practical OPE is insufficient to reproduce the experimentally fa- 
vored shapes of DAs |200| . The same conclusion can be reached by a comparison of a re-ordering of the 
OPE, which sums up an infinite series of operators of increasing powers in covariant derivatives, with 
experimental data. 



7.1 Distribution amplitudes and photon-photon annihilation 



Before going into the details of the sum-rule approaches of |185j and |2()()j to the determination of DAs 
we will set the stage by explaining and applying them in (3-1-1) dimensions. In the 2D case we have 
already encountered them in our discussion of the 't Hooft model. 

In the calculation of exclusive processes with a large momentum transfer - like the hadronic decay 
of heavy mesons or the photon-photon annihilation into hadrons - mesonic DAs appear quite naturally. 
They were introduced to separate nonperturbative large-distance physics from the perturbatively ac- 
cessible small-distance physics in the calculation of transition amplitudes 11901 HE?] (an exhaustive 
discussion is given in jl9Up . In general, the DA of a hadron decomposes into a part describing the valence 
quark and/or gluon content and into parts associated with the contribution of higher, color-singlet Fock 
states. For example, the 7r+ wave function ip{x,fi'^) associated with the lowest-twist, nonlocal operator 
and with the lowest Fock state is defined as 



(i(z)7^75Pexp 

n=0 "'■ 

iq^ (j){z ■ q, fi^) + ■ ■ ■ , 



ui-z 



TX 



dizh,^, [z'^i'S^ - S^])"m(0)| 7r+(g) 



(76) 



where the separation 2z is light-like, z^ = 0, the sum is over even n, and V, fi denote path ordering 
and the normalization point of the operators, respectively. In Eq. fl76|) only the function (p^iz ■ qjfJ-"^) 
associated with the valence quark distribution is explicitly indicated, paranthesis refer to higher Fock 
states. Writing ^ = 2x — 1, the quantity x or {x — 1) is the fraction of longitudinal pion momentum 
carried by the quark or antiquark, respectively. The DA 99^(^,/i^) is interpreted as the integral of the 
full Bethe-Salpeter wave function over transverse momenta k± down to k± ~ /i. The dependence on the 
cutoff II in (fni^, yU^) can most conveniently be obtained by expanding into eigenf unctions of the evolution 
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Figure 16: Diagram contributing to the hard scattering amphtude in 77* — Tr*^ to lowest 
order in as- Wiggly hues denote photons, a sohd hne a quark or an antiquark. 



kernel obtained by considering one-gluon-exchange ladder diagrams which contribute perturbatively to 
the wave function |186j . The expansion is given as 



n=0 



log 



A2 



-7n 



(77) 



where the eigenvalue 7„ is given as 7^ = 3^ 1 — 2/{{n + l){n + 2)) +4X]r=2 1/^ ■> ^"^^ ^ = -'-^ ^ 3 
is the first coefficient of the QCD beta function. Due to the ort honor mality (with weight (1 — ^'^)) 
of the Gegenbauer polynomials C^^^ on [-1,1] the expansion coefficients are simply given as a„ = 
j\d^ C^/^(^)(y9(^, /Iq). Let us now discuss two applications of DAs. We first consider the annihilation 
of two photons into hadrons, 77 — ^ Xh- One of the photons either is highly virtual or it is real but 
incident under a large angle Oc.m. in the center-of-mass frame. This is described in terms of a convolution 
of a hard scattering amplitude Th, calculable in perturbation theory (see Figs. ll6ir7j) . with the DAs 



of the hadrons produced. For example, for the process 77* 



this arises by considering a partial. 



perturbative contraction of quark fields in the amplitude (7r°|j^(?/)j^(0)|0), yielding propagator lines to 
which the external photons couple, and a soft remainder which in lowest-order twist is given by Eq. (|77|). 
The vertex for this process is 



(7J 



where p is the pion's momentum, the polarization vector of the real photon, q (with = — Q^) the 
momentum of the virtual photon, and F^^{Q^) the photon-pion transition form-factor. The latter can 
be written as the convolution of hard scattering amplitude (see Fig.lTBj) with the pion wave function. 
To zeroth order in ag we have 



xil — X 



(79) 



where = min(a;, 1 — x) Q. The amplitude A^aa' for the process 77 — > vr+vr at large angle 9c.m. is 
given as |189j 



A^AA' = Jq ^^2¥''r(a;i,(5^Jv^,r(a;2,<5^2) x Txy{xi,X2, s,9c.m.) 



(80) 



where A, A' denote the helicities + or — of the photons and Q^-^ ~ min(xi, 1 — Xl)^/s\ sin^^.m. I (accord- 
ingly for Qxi)- III the case of equal helicities one has = M. which can be expressed in terms of 

^■^ Abusing the notation, wc denote the different functional dependences of the wave function on 1^ and x by the same 
symbol ipT^. 
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Figure 17: Diagrams contributing to Tw of Eq. (jSUI) to lowest order in a^. Wiggly lines 
denote photons, a solid line a quark or an antiquarks, and a dashed line is a gluon. 



the pion's electromagnetic form factor as 

1 - COS^ 6c.m. 

where 



3s Jo JO xi(l — Xi)x2(l — a;2 

Notice that the expansion of the form factor F^^ in Eq. (|82p starts at linear order due to the kinematic 
situation shown in Fig. El - in contrast to the process 77* 7r° (see again Fig.ll6|). 

If the respective form factors are measured accurately - like the CLEO collaboration ' 191j did in the 
case of -FV^(Q^) - then Eqs. (fTUIl . ljH^ can be used to constrain the form of the pion DA. For Q — 00 it 
should approach the asymptotic parton-model form V5^(x, Q^. — > 00) = 6/^x(l — x). 



7.2 OPE and OAs 

7.2.1 Pion DA and local condensates 

Armed with the concept of DAs and the conviction of its usefulness we are now in a position to relate 
them to the OPE parameters characterizing the nonperturbative QCD vacuum - the gauge invariant 
condensates. This idea was pioneered by Chernyak and Zhitnitsky |184[ 1185] (for a review see [1 82| ) • 
We will mostly follow their presentation which focuses in particular on pion wave functions. 

As one can derive from Eq. (j76|) . the following relation holds for the moments of the wave function, 

(e% = /-\'c^er(^(e,/i^), 

^fA^ ■ qf^' (e)^, = (0|d(0)7.^'^75 {^z^0 , - ^D,])"" «(0)|7r+(g))^^ (83) 

where (p is a dimensionless wave function, (f = j-ip. According to Eq. ()83|) we have for the zeroth 
moment: 

{e) = 1 . (84) 

To make contact with vacuum condensates the central object to consider is the correlator 

= {z-qf^^'^^^'lr^Mq') (85) 
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of the local operators 



On = ^u-f^z-'-f^ (iz'^l'S^ - 'D^]yu- (u^d). (86) 



Expanding IniN2{.t) iiito an OPE in the euclidean region q = —Q < 0, we have 



(Ar + l)(Ar + 3)47r2 ^V/^V 12(g2)2\7r 

(l + ^ [(iVi - N,f -N]) + {uu)l, (ll + 4iV) (87) 



1 



where N = Ni + N2, and SU(3) p symmetry and vacuum saturation for the four-quark condensates have 
been assumed. The index at the dimension-six condensate refers to the normalization point. The 
intermediate states in the correlators Iqq and Inq are the same, the resonances to which the "currents" 
Oq and Otv both couple are the 7r+ and spin-one excitations such as the Ai, ■ ■ ■. We assume the first 
to be explicit, and we suppose that the rest can be treated as a part of the perturbative, spectral 
continuum which starts at some threshold sno- In the chiral limit m^r — we obtain the following 
Borel sum rule for the correlator 

f/^^\ = ^ (i ^ e-'No/M^] _^ 

^-X*^ /m2 47r2(Ar + l)(Ar + 3) V 



12M2 + 1 \ TT / 81M4 

Evaluating the sum rule in Eq. ()88p for the first three moments = 2, 4, 6 by using the SVZ values for 
the condensatesl^ (^G^) = 1.2 x IQ-^ GeV^ and a, (uu) = 1.83 x 10"^ GeV^ the following resuhs at 
fil = (0.5 GeV)2 were obtained in |185j 



n.... X. =0-46, ., = 0.30, ..=0.21. (89) 



= 46 /f^^ 

/(0.5GeV)2 ■ ' \^ /(0.5GoV)2 ' \'' /(0.5GeV)2 

Using the convex asymptotic form (fas = |(1 ~ ^^) for /i^ — > cc - the contribution at n = in Eq. (f77j) 
- one arrives at too low moments in comparison with Eq. ()89|) : = 0.20 , (^^)^^ = 0.086, and 

{^^)as = 0.048 . This form of the wave function seems to be excluded by the sum rules for the moments. 
Therefore, Chernyak and Zhitnitsky proposed the following pion DA 

^cz(e, f^l = (0.5 GeV)2) = ^ - e) (90) 

which is not convex anymore, strongly concentrated at the endpoints ^ = ±1. This DA gives moments 
that are in reasonable agreement with Eqs. 



f)^^ =1, =0.43, U'^)^^ =0.24, U^)^^ =0.15. (91) 

^ /(0.5GcV)2 ' \^ /(0.5GoV)2 ' \^ /(0.5GoV)2 ' \^ /(0.5GoV)2 ^ ^ 

Working with ipcz to describe branching ratios of pionic decays of charmonium levels, experimentally 
obtained values can well be matched |185j . In addition, it can be demonstrated using i^cz that the 
annihilation contributions to mesonic D decays are comparable to the direct ones leading to a better 
agreement between theory and experiment. 



Anomalous dimensions in the four-quark operators have been neglected, for a discussion of the evaluation of the sum 
rules see | il85j . the SVZ value for the condensates it taken at /x = 1 GeV. 
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7.2.2 Pion DA and nonlocal condensates 



The approach of Chernyak and Zhitnitsky, which relates the moments of DAs to the local condensates 
appearing in the OPE of a suitable current correlator (see Eq. (jHS))), implicitly assumes that quarks and 
gluons in the nonperturbative QCD vacuum have zero momentum. If this was the case for quarks, the 



nonlocality parameter A'^, defined as 



A^ 

q 



(g(O)g(O)) 



(92) 



q denoting a light-quark field, would vanish QCD sum rules, however, predict in the chiral limit a 
value of Ag ~ 0.4 — 0.5 GeV^ with an error of 10-20% j l93tll94j . An estimate using the single-instanton 
approximation of the instanton liquid yields the slightly higher value A^ ~ 0.6 GeV^ |195| il96j while an 
unquenched lattice computation with four fermion flavors yields values close to the sum-rule estimate 
|248| I244j . A more recent investigation of lattice estimates and the extrapolation to the chiral limit 
thereof was performed in |182j . As a result, a window 0.37 GeV^ < A^ < 0.55 GeV^ was obtained. 
vPhenomenologically, it is thus an established fact that gauge invariant correlations between quark 
fields are of finite range. Moreover, the mass scale governing the fall-off of the quark correlator in 
Euclidean spacetime is larger than the usually accepted condensate scale ~ Aqcd ~ 0.3 — 0.4 GeV which 
makes the OPE philosophy of an expansion in powers of the ratio Kqcd/Q questionable. In relation 
to the moments of DAs the problem can be quantified as follows. According to Eq. (j88p the ratio of 
nonperturbative contributions to the perturbative part in the OPE is (at large A^) cubically growing with 
A^. This means that even in a truncated OPE, effectively, the nonperturbative scales in the vacuum, 
which determine the A^th moment of the DA, appear to be rapidly growing with A^. To compensate 
for this effect higher and higher values of the continuum threshold Sqtv have to be adopted to reach 
stability of the sum rule for the A^th moment |185j . But this amounts to neglecting higher resonances 
in the spectrum that would be related to the OPE contributions governed by large, nonperturbative 
mass scales. The quality of such a cancellation between neglected parts of the integrated spectrum and 
neglected condensates of a large mass scale in the OPE is obscure. It is very reasonable though to 
expect it to become worse with growing A^. Considerations of this type lead Mikhailov, Radyushkin, 
and Bakulev to postulate a modification of the OPE with built-in nonlocal effects |197t I19H1 11991 I^UU] . 
An approach, which expands nonlocal condensates into complete sets of resolution dependent, finite- 
width functions in position space, was proposed in and will be discussed in detail in Sec. 17.3.21 with 
respect to the question of local duality violation. Instead of writing sum rules for the moments {^^^^ 
as in Eq. (|HH|) one can use the truncated OPE to directly write a sum rule for the wave function ipT,{x) 



4^ 



(1 



7r«, {qqf (11 [6{x) + 5(1 - x)] + 2 [6' {x) + 6' {I - x)]) 



81M4 



(93) 



where the prime denotes differentiation with respect to x. Taking into account radiative corrections in 
order a^, it was obtained in ^200j that the foUowinf substitution should be performed 



iPn{x,fJ.^ 



oo 



) (Pt,{x,ix^ oo) 




(94) 



^^The parameterA^ is a measure for the deviation from perfect, gauge invariant field correlation at different points in 
spacetime. It appears in the first nontrivial term when expanding the corresponding gauge invariant two-point function 
(g(z)7'exp \ig dy^ ^niu)] ^(O)) in powers of z^. 
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Note that this form of the DA, which arises from a truncated local expansion of the correlator of 
nonlocal, gauge invariant currents, violates the boundary condition 

'^{x,iJ,^) < Kx' for X ^ (e > 0) , (95) 

which assures the convergence of the expansion in Eq. ()77|) |186j . Hence Eq. ()94|) is incomplete. The 
complete form of Eq. (jHSD is obtained by substituting 5^5 + 21og(MV^2) ^q. (jHU piHI^ . 

The starting point in |197| 11981 1199| I2UU1 l2U5j is a parametrization of the gauge invariant bilinear 
scalar or vector quark condensate (or correlator) of the form 

(g(O)g(z)) = / e^-'/'fsiiy) , (g(0)7^g(^)) = zz^-na,{M') {qqf / e'^^'/'Mu) (96) 
JO 81 JO 

where fs and fv denote the distribution function for quark momenta in the vacuum associated to the 
scalar and vector correlation on the left-hand sides, respectively. In writing Eq. ()96|) the Fock-Schwinger 
gauge together with a straight-line connection between the point and z were assumed. Higher n-point 
functions as well as the nonlocal gluon condensate are parametrized in a similar manner, see \2[)(}\ l2Ulj 
for an extended discussion 

The functions fs and fv have to be modelled. For condensates without covariant derivatives in the 
OPE a formal expansion of the functions fs, fv into the set 5^"'\u) , (n = 0, 1 ■ ■ ■) is truncated at = 0. 
Operators in the OPE with n powers of covariant derivatives are associated with the 5"(i/) part of the 
expansion. We have already seen that there are problems with this local truncation. Mikhailov and 
Radyushkin proposed to model fs and fv by simply shifting the argument of the delta function of the 
local expansion to generate finite-width correlations in position space 

^ _ ^l^) . (97) 

In the case of higher n-point functions the / functions are modelled as products of delta functions 
centered at nonzero values of z/j. As a consequence, the expression for the pionic DA is smooth at x = 
and the moments of the nonperturbative terms decrease with in such a way that the ratio to the 
perturbative terms increases only moderately - much in contrast to the case of a truncated OPE, see 
Eq. (I88|) . Other ansatze than the one in Eq. (|97|). which are better suited for sum rules with nondiagonal 
correlators, were proposed and discussed in |2(J2| 12(13] I2(J4| ITH^ . Using A^ = 0.4 GeV^, the following 
values for the first three moments were obtained 

(^2) = 0.25 , (^^) = 0.12 , (^^) = 0.07 , (98) 

which are close to the asymptotic values, see text below Eq. ()89|) . For recently obtained small alterations 
of these values see |223j . These values can be well reproduced by using the model DA = ^^x{l — x) 
which differs quite drastically from the one in Eq. ()90|) obtained by Chernyak and Zhitnitsky and is 
qualitatively similar to the asymptotic one. 

7.2.3 Direct experimental evidence for nonlocal condensates from meson DAs 

The only way we can decide whether at a given order in twist and the prediction of an endpoint 
concentrated pionic DA, which rest on a truncated OPE (Chernyak and Zhitnitsky), or a prediction 
close to the asymptotic, convex form, which is based on nonlocal condensates (Mikhailov, Radyushkin, 
and Bakulev), is true is to compare them with precise and independent experimental data on exclusive 

^^In the case of the gluon correlator <^ylJ^(z)A|j,(y)) (with straight lines connecting the origin with z and y) one would 
in principle not only have a dependence on (z — y)'^ but also on and y^. However, a local expansion reveals that the 
dependences on and y^ are much weaker than the one on [z — y)^, and therefore one may neglect them. 
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quantities. A very well suited observable is the pion-photon transition form factor F-,j^[Q'^) of Eq. ()79|) . 
It was measured rather precisely by the CLEO collaboration for from 1.5 to 9 GeV^ |191j . Together 
with the data taken by the CELLO collaboration |192j a range of values from 0.7 to 9 GeV^ is 
covered. 

In |206J a perturbative approach was applied to the calculation of the transition form factor Ft^^{Q'^) 
which takes into account transverse momentum effects and the Sudakov suppression of longitudinal 
momenta close to the endpoints x = 0, 1. More precisely, the (perturbatively undressed) valence Fock 
state pion wave function \E'o(x, is assumed to factorize into the DA ip{x,fi'^) for the longitudinal 



momentum fraction and a portion — x)b) containing information about the distribution in two 

dimensional transverse position space 



vI/o(x, 6, /i^) = ^ <^(x, /i^) S( - x)h) (99) 

where /x is the scale for the factorization of hard scattering and soft momentum distribution, and the 2d 
vector h denotes the quark-antiquark separation. Certain constraints on E can be derived from duality 
arguments |210j which are minimally satisfied by the following Gaussian ansatz 



— x)h) = Att exp 



x(l — x)b'^ 
4^2 



(100) 



where a denotes the transverse size parameter. It can be fixed by the requirement that the wave function 
'^o{x,b, fx'^) be normalized to y/G/fn (due to 7r° — >• 77). For models with a power-law ansatz for the 
dependence of the wave function on transverse momentum see |21H I212[ I213| I214j . Taking into account 
Sudakov suppressions in the hard scattering amplitude Th when retaining its dependence on transverse 
momentum k± |215^ \'2W\ I217j . the expression for Ft^^^Q"^) is altered in comparison to Eq. (ffnjl as 

F^^iQ^) = C dx ^^o{x, b, = 1/6^) Th{x, b, Q^) exp[-S{x, b, Q^)] (101) 
JO 47r 

where S denotes the Sudakov exponent which accounts for gluonic radiative corrections not contained 
in the evolution of the DA (f. Note that the factor in Eq. (jl01|) represents an ansatz designed to 
take effects of intrinsic transverse momentum into account, it is not derived from first principles. The 
accuracy of this ansatz was questioned in |207j . Also, subleading logarithms in the Sudakov corrections 
were taken into account in |207j . Still, using the asymptotic form = 6x(l — x), very good agreement 
between the data and the numerical result based on Eq. pOlj) is obtained at /x^ = 1 GeV^ in ^^i^, see 
Fig. (|18|). A perturbative investigation of the transition form factor 

FMQ^Q'") = % f'dx /[""^ (102) 
V3 "'o xQ-' + (1 - x)Q' 

for the process 7*7* 7r° using the same approach as in |206j was performed in |218j . At ~ Q'^ this 
calculation provides at test of the hard scattering amplitude since the transition form factor is essentially 
independent of the pionic DA. The result of |206j that the DA of Chernyak and Zhitnitsky is excluded 
and the asymptotic DA is favored by experimental data was confirmed in |218j when taking the limit 
Q'^ 0. However, a more recent light-cone sum-rule analysis of the CLEO data |208j indicates that 
the neither the asymptotic DA (although favored) nor the Chernyak- Zhitnitsky model are compatible 
with the data. This was qualitatively confirmed in |209j where an incompatibility with the CLEO data 
on the 3 a and 4 a level was obtained when using the asymptotic DA and the Chernyak-Zhitnitsky 
model, respectively. In |219j the transition form factor F^^* (Q^, Q'^) was calculated using a QCD sum 
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The analysis in |2U8| obtains constraints on the first two Gegenbauer polynomials. 
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Figure 18: Comparison of the CLEO-results (points) for Q'^F-^^^^Q^) with the theoretical pre- 
dictions made by |2U6j using the asymptotic DA (sohd curve) and the Chernyak-Zhitnitsky 
DA of Eq. (j9(Jj) (dashed curve). The dotted curve shows the prediction made with the 
Chernyak-Zhitnitsky DA when its QCD evolution is taken into account. The solid line 
at 2/^ indicates the asymptotic value. Plot taken from |191j . 



rule for the three-point correlator associated with the process 7*7* vr°, see Fig. for the lowest- 
order perturbative diagrams contributing to it. Applying a factorization procedure to separate off the 
infrared singularities emerging in the limit Q'^ — 0, subsequently absorbing these singularities into 
bilocal correlation functions, and finally assuming certain forms for these correlation functions, the 
dependence of the transition form factor Ft^^{Q'^) was studied and compared with the data from CELLO 
|192j . see Fig.|2ni As a result, again evidence was provided by this approach to F.,,y{Q^) that (f{x,Q^) 
is rather close to the asymptotic form (f{x, — >■ 00) = 6x{l — x). The same conclusion was reached 
in [220) where a light-cone sum rule (see next section for its foundations) for the two-point correlator 

d'xe-''^' {7c\p)\Tj,{x)j.{0)\0) (103) 

of the electromagnetic currents was used in the variable^® Q'^. In the chiral limit the light-cone OPE 
for the transition form factor F^^^IQ"^) can then be written in terms of a dispersion integral over the 
imaginary part of the transition form factor Im F^^^y* {Q'^ , s) which, in turn, is given by the pion DA up 
to twist four as 



-\mF^^*[Q ,s) - 



TT -7 V- ' 3 I s + Q2 Q2 



(104) 



where ^p^^^ denotes the twist-four part. As a result, even at a higher order in twist the asymptotic DA 
is again favored by the CLEO and CELLO data, see Fig.|2Il Notice, however, that a very recently 



^^Photon virtualities are again — and Q'^ = — (p — 



41 




Figure 19: Lowest-order diagrams in perturbative QCD for the three-point correlator 
27ii J d'^xd^y e~^'^^'e^P^ {Tjl{y)j^{x)ju{0)) where jf = U'y^'jxu — d'-y^'yxd, and denotes the 
electromagnetic current, = |u7^m — ^djf_id. Taken from j219j . 



performed analysis of the CLEO and CELLO data using a light-cone sum rule with the inclusion of 
twist four and 0{as) perturbative corrections to the pion DA seems to exclude both the asymptotic as 
well as the Chernyak-Zhitnitsky form of the DA at lowest twist \2'2'2\ 12211 . Still, in the determination 
of the pion DA from a nonlocal sum rule at the same order in twist and as a nonlocality parameter 
Aq = 0.4 GeV^ is compatible with the result obtained from the light-cone sum rule. Seemingly, the 
inclusion of radiative corrections has a greater effect on the deviation from the asjTuptotic DA than 
the inclusion of the twist-four correction. This fits well with the result of |223j where the sum rule for 
the moments (^^^'^ with nonlocal condensates (including the Ai resonance explicitly into the spectral 
function) and 0{as) corrections but no twist-four corrections was analyzed. A number of calculations of 
the twist-two pion DA (f were recently performed using instanton models |225| I226j or in an instanton 
motivated model |227j . As compared to the asymptotic DA all of these investigations found a slight 
depletion of the model DA within the central region but no qualitative deviation from the asymptotic 
form. 

We may thus conclude that the compatibility (noncompatibility) of the pion DA directly extracted 
from the experimentally well measured transition form factor F^^jlQ'^) and the pion DA obtained from 
a sum rule involving a truncated, practical OPE with nonlocal (local) condensates provides rather 
strong phenomenological evidence that a local expansion of nonperturbative gauge invariant correlation 
functions is not sufficient to describe the hard and exclusive process 77* —>■ ir^. 

The technically more involved and theoretically more shaky case of p meson DAs will not be discussed 
in detail in this review^^. We only mention here the QCD sum-rule analysis of Refs. |229| l2U5j using a 
truncated, practical OPE with local and nonlocal condensates (A^ ~ 0.4 GeV^), respectively. In both 
cases only leading twist was considered and, in contrast to the old calculation of Refs. |185| 1233] I234j . 
radiative corrections to the perturbative parts of the sum rules were taken into account. In |229j a 
sign error in the four-quark contribution to the OPE for the second moment of was pointed out 
leading to a drastic change in the prediction of the shape of this DA as compared to the prediction 
in \1E51 \23'3. I^4|j . A comparison between (p\\ at a normalization scale fi = 1 GeV, obtained from the 
sum rules for the zero-helicity state reveals no essential difference in both approaches. The treatment 
with nonlocal condensates, however, gives a much better stability of the Borel curves of the moments 
(^^) up to large N than it does in the local case where a large- prediction is unreliable, see Eq. ()88p 
for a similar situation in the pion case. The shape of (f^ obtained in |229t I205j is rather close to the 

^^Instead of one DA there are four DAs (p±,if^^,g^, and at the leading twist two |228[ 1229] and valence Fock state. 
The DAs ip± and describe the distribution of longitudinal quark momentum in a transversely and longitudinally 
polarized p meson, respectively. The twist-two DAs and 5° are connected to (^|| |229j by Wandzura-Wilczek type 
relations (230) whose status as a dynamical statement is still under debate, see for example [2321 1231j . 



42 



0.30 



0970597-011 



□ CELLO 
♦ CLEO 



2.5 5.0 7.5 10.0 

(GeV^) 

Figure 20: Comparison of the CLEO-results (points) for Q'^F^.^^Q^) with the prediction of 
[219 which is based on a QCD sum rule for F^^. (Q^ = —q'^,Q'^ = —q'"^) with two virtual 
photons. The limit Q'^ was taken after applying a factorization procedure for the 
infrared singularities which arise perturbatively in this limit. Taken from |191j . 



result in |185t 12331 1234j which is not too far from the asymptotic one; for <^_l, however |229j get a much 
wider distribution than Chernyak and Zhitnitsky. Let us also mention that a computation of ip± was 
performed in the framework QCD sum rules with nonlocal condensates in |224j . 

To summarize, we have seen that the use of nonlocal condensates in the OPEs of current correlators 
which determine the moments of mesonic DAs is superior to the treatment with local condensates in 
the following ways: (i) the prediction of high moments (^^^^ is feasible, (ii) Borel sum rules for a given 
moment exhibit a wide window of stability, (iii) the prediction of close-to-asymptotic behavior of the 
lowest-twist pion DA using nonlocal condensates in spectral sum rules for the first few moments agrees 
well with an independent determination in terms of light-cone sum rules using the experimental data 
on the transition form factor F^^(Q^) as a phenomenological input. 

7.3 Reshujflmg the OPE 

Motivated by the phenomenological results in the last section we will in this section review the theoretical 
foundations for a light-cone expansion of a current-current correlator into gauge invariant, nonlocal 
operators. Taking the hadron-to-hadron (hadron-to-vacuum) matrix elements of this expansion allows 
to relate experimentally measurable deep inelastic scattering (transition) cross sections to structure 
functions (distribution amphtudes). 

A phenomenological expansion of (2-point) nonlocal operator averages and Wilson coefficients into 
nonlocal objects, which captures more information about the nonlocal operator average than its trun- 
cated local expansion, is proposed for the vacuum-to-vacuum case. This approach introduces a so-called 
resolution parameter. The nonperturbative evolution of the nonlocal objects in this parameter will be 
derived. Several applications will be discussed, and a "running" gluon condensate will be extracted 
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Figure 21: Comparison of Q"^ F-n^^lQ"^) as obtained from the analysis using a light-cone sum rule |22()j 
and the experimental data. The solid line is associated with the asymptotic DA, the long-dashed 
line with the Chernyak-Zhitnitsky DA, the dash-dotted line with the simple interpolation formula 



and the short-dashed line with a DA obtained in 



Plot taken from 



from experimental data in two channels. The implications of our nonlocal modification of the OPE in 
view of the nature of this expansion and an OPE-based realization of local quark-hadron duality are 
investigated phenomenologically. 



7.3.1 Light-cone expansion into string operators of a current-current product 



We have already mentioned in the last section the use of light-cone sum rules as a means to extract 
information on the twist-two pion DA from the experimental data for the transition form factor 
|220j . In this section we will give some theoretical background for this expansion. 

The light-cone expansion of an object like iTj^{x)ju{—x), jfj, being the electromagnetic current, into 
nonlocal string operators |24()| 1239^ l2Hf)j and the know-how about the perturbative renormalization of 
the latter under a change of scale are of paramount importance for the experimental determination of 
hadronic structure functions (hadron-to-hadron matrix element of nonlocal operator) and DAs (hadron- 
to- vacuum matric element of nonlocal operator). Within perturbation theory the proof for the existence 
of the light-cone expansion of a current-current correlator in terms of quark twist-two string operators 
hke^o 



'ijj{x)X°'xV exp 



^(0) =^{x)X'x[x,Q\i){Q) 



(105) 



was given in |235j . In Eq. ()105p A° is a flavor Gell-Mann matrix and x = '^^x^. Some time later the 
one-loop dependence of gauge-string operators on the renormalization scale /z was calculated in j 236j by 
integrating over fluctuations q'^ and A'j^ about classical backgrounds q'^ and A^. The relevant diagrams 
are shown in Fig. 1221 

Using the Fock-Schwinger gauge x^A^^ = 0, the background gauge-string [ax, 0] is unity and thus 
is omitted in some of the expressions to follow below. Evaluating the diagrams in Fig.|221 and keeping 
only the first logarithmic term in x^, one then obtains (236|| 



ip{x)X''xip{0) 



da / d[3 




5{l-a)5(3 



^log4i^(«,/5) 
Stt /if 



ip{ax)\°'xip{l3x) 



(106) 



^"^We consider the flavor nonsinglet string operator first since for this case we do not have to consider the mixing with 
the gluon string operator. 
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Figure 22: Diagrams contributing to the one- loop renormalization of the flavor- no nsinglet 
string operator of Eq. (jlOSj) . Thick lines denote the classical gauge string, crossed lines 
classical fermion fields and dashed lines are fluctuating gluon fields. 



where 



1-/3 



-6{P) / dp' 



P' 



-m 



a 



1 — a 



— 5{a) / da 



a 



1-a' 



The leading logarithms in Eq. 



(107) 

can be summed by solving a one-loop renormalization-group equa- 
tion. In analogy to an expansion of a local operator into an eigenbasis (conformal operators, multi- 
plicatively renormalized) of the evolution kernel one may also expand a nonlocal string operator into a 
nonlocal eigenbasis. The result is[ 



a + P a — P 
X, X 



' a — P 



dk 



-ika 



X ifj 



Air 

T fa' + p' 



l/2+ioo 1 

dju + -)p-'/'j,+i/2m 

— 1/2— ioo Z 



Ma 

oo 



dp'^P'H]^,^,{kp') 



X A'^x 



a' + P' a'-P' 
X, X 



'a'-P' 



-X 



where Jj+i/2(-f^j+i/2) is a Bessel(Hankel) function. 



7i 



36 Jo 



f^duu^'^ -6(l-u)-(l-u)-2(— 5(u) I du'^^ 

Jo 2 \ 1 — u J 1 — u 



'a(jj^ 



-7j 



:io8) 



(109) 



and the integrations over a', P' and k,j indicate the projection onto and the (continuous) decomposition 
into renormalization-group covariant, conformal string operators, respectively. 

In addition to the quark string operators of Eq. (jlOSp there is a twist-two flavor singlet gluon string 
operator which enters the light-cone expansion of a current-current correlator. It is given as 



G{U,V) = XaG"^^{ux) [uX,Vx\abG^^p{vx)xp 



(110) 



where the straight-path Wilson line [ux,v,x]ab is now in the adjoint representation! The flavor singlet 
quark string operator 



Q{u,v) 



ip{ux)x[ux, vx\iIj{vx) — iIj{vx)x[vx, ux\iij{ux) 



;iii) 
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Figure 23: Diagrams responsible for the mixing of twist-two, flavor singlet quark string 
operators and flavor singlet gluon string operators. 

and the gluon string operator of Eq. (jllUj) mix under renormalization, see Fig. (j23|) . Formulas simular 
to the one in Eq. (jlU8|) relate Q and G at a renormalization point to a linear combination of their 
respective conformal expansions at some other renormalization point /if, details can be found in |236j . 
The derivation for the evolution of string operators containing three Wilson lines is similar except for 
the fact that conformal symmetry is not sufficient to determine the solution explicitly. 

It was shown in |23tit I239j how the light-cone expansion of a T product like Tj^{x)j^{—x) into 
nonlocal string operators instead of local operators (in the latter case the choice of a suitable operator 
bases at a given nonleading twist is obscure because there are relations between the local operators 
dictated by the equations of motion |238t 1237] ) leads to a compact and much better managable series 
in singularities in deviations from the light-cone = than in the local case. The leading term 
— -^^^2^4 V'(^)7m^7i^Q^^(~3:) , where Q = + ^^A^ is the charge matrix, has no definite twist. To 
project onto the leading twist, one can perform a symmetrization and trace subtraction in the local 
expansion of this string operator, which can be reassembled into a nonlocal expression. Higher twist 
contributions are technically quite involved. 

The nonlocal light-cone expansion of |236j was used extensively in phenomenology. Besides the 
applications already addressed in the last section let us just mention two more examples: Light-cone 
sum rules were employed in |241j to predict the heavy-mesons couplings to pions in nonleptonic D and 
B decays in terms of two- and three- particle DAs up to twist four. The pion form factor at intermediate 
momentum Q ~ 1 GeV was estimated in |242j using light-cone sum rules with higher twist distributions 
and radiative corrections. 

7.3.2 Delocalized operator expansion 

After having discussed the expansion of a T product of electromagnetic currents into nonlocal string 
operators in view of applications involving a nonlocal hadron-to-vacuum matrix element or a nonlocal 
hadron-to-hadron matrix element we turn now to the case of a vacuum-to-vacuum average which leads 
to the occurrence of nonlocal condensates. The discussion of [7^], which we follow, will be on a more 
phenomenological level. We have already introduced the approach of |197| 11981 ll99U2U0ll205j to nonlocal 
condensates in Euclidean spacetime in Sec. 17.2.21 

Let us first discuss the case of two-point condensates. They are associated with nonlocal versions of 
the quark and gluon condesates in the conventional OPE. The parametrization of the nonlocal scalar 
quark condensate in Eq. (pUj) is to a good approximation also applicable to the nonlocal gluon condensate 
in Euclidean spacetime: In Fock-Schwinger gauge we have 120 Oj 

{A;{z)K{y)) = 6'^' {y,z^ - 5,^{z ■ 1/)) ^ x g{{z - y)\ z\ y') + ■ ■ ■ (112) 
where the local expansion of g{{z — y)^, z^, y"^) reads 

,((. - y)\ z\ y-) = l- ^'"^ ((y - zf + + • ■ ■ - ,((. - yf) . (113) 



46 




Figure 24: Schematic drawing of the short-distance function f{x) and the 2-point condensate 
g{x) illustrating the scale hierarchy A//Ag~A/(5-Cl. Plot taken from [73] . 



In Eq. ()112j) the paranthesis denote contributions to (^A'^{z)A''j^{y)^ which are not captured by the scalar 
function g{{z — yY, z^^y"^). Clearly, the dependence on {y- 



, y^). Clearly, the dependence on {y — zY is much stronger than the dependence on 
y^ or z^. The dependence on the latter two variables thus can be omitted on the level of accuracy at which 
the local condensates are known. In the expansion of the scalar part T(g^) of a given current-current 
correlator i J d^x (Tj^(a;)j'j^(0)) the Lorentz, color, and flavor indices in front of the nonlocal condensate 
get contracted and integrations over x, ?/, and z are performed. The nonperturbative correction to 
Tperti^f) arising from a two-point condensate g{x'^) can thus be written as 



d'^x f{xi,. . . ,Xi) gix"^) 



fni,...,n4,i^) 9ni,...,n4, 

ni,...,n4=0 



(114) 



where 



fni,...,n4,{^) 



d^x f{xi, ...,Xi) 



n 



d'^x 
d^k 



d 



4 Qrii+l 

ll^^H^X^x, 

i=l 



n 

i=l 



nil 



e 



(115) 



fci,...,A;4=0 



~9{k 



g{xi, . . . ,X4) 



(116) 



A number of comments are in order. In Eq. (jll4p the function /(xi, . . . , X4) denotes the strongly-peaked 
at a; = 0, perturbatively calculable, short-distance dependence as it arises in Euclidean spacetime after 
separating off the condensate part in an application of the background- field method^^. It corresponds to 
the Wilson coefficient for the local condensate in position space and has a width A f which is comparable 
to the inverse external momentum Q~^, see Fig. (j2l|l . To arrive at the second line a complete set of 



'^We have notationally suppressed the dependence of / on the external momentum Q. 
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functions e^^...„^(x) and e^^...„^(x), which span a dual space with respect to the bihnear form 

if,9) = j d'^xf{xi,---,Xi)g{x^) (117) 

and which are normalized to 

( ^ni,---,n4 5 ^mi,---,m,4 ) ~ ^nimi ' ' ' ^runn i (H^) 

has been inserted inbetween / and g. In the hmit, where the resolution parameter Q goes to infinity, 
Q ^ oo, the set e^j^...n^{x) and e^^...„4(a;) reproduces the local expansion of the nonlocal condensate, 
and we have 



V 5("^)(xi)---5("*)(x4) and e^;:.Zi^) = x'l' ■ ■ ■ ■ (119) 



...'n ,\X j 



^l! ■ ■ ■ n^l 

In Eq. pi9j) 5*^"') denotes the n^th derivative of the delta function. For definiteness, we have chosen in 
Eqs. fjllSp . (11161) the following set of products of Hermite polynomials and their duals 

Q{ni+l)+-+(n4+l) 



Obviously, the basis ()120p reduces to the basis pi9p in the limit Q —>■ oo. Different choices of one- 
parameter basis functions are equally well possible. Let us emphasize at this point that the delocalization 
of the OPE, which is achieved in this manner differs from |197| I198| 11991 I^UUt l2U5j by the fact that 
with a truncation of the expansion (jll4|) at ni = 0, ■ ■ ■ , n4 = a resolution parameter, which eventually 
will be chosen to be equal to the applied external momentum, cuts off the irrelevant long-distance 
contribution of the nonlocal condensate gix"^) to the integral. The Wilson coefficient in momentum 
space is equal to the one of the local expansion in this truncation, and this leads effectively to the 
running of the condensate with resolution. In Sec. 17. 3. 71 an alternative approach to running condensates 
will be discussed and applied. 

The properties of Hermite polynomials dictate the following evolution equations for /ni,...,n4(^) and 

Af (O] - KulK. , I K-l)n4 

Jni,...,n4\^^J — 2il^ Jni-2,...,r!,4l^'J "T • • • "T 2 J ni,...,n4-2\^ ^ J , 

An (0\ - (ni + l)(ni+2) ^ 

9ni,...,ni{^^) — 20,^ yrJi+2,...,n4l^'J 

K + l)(n4 + 2) 

■■• 5'ni,...,n4+2li'J • l-L^J-J 

The parametric counting of the short-distance coefficients /ni,...,n4 and the long-distance functions 
gni,...,n4 is similar as in the local expansion for Q > Q, the leading local terms are corrected by a 
finite sum over powers in Q/Q in the former case and an infinite sum over powers of A/f2 in the latter 
case. Here A denotes the mass scale associated with the fall-off of g{x'^) in Euclidean spacetime, for 
details see Sec. fj7.3.3p . 

The consideration of 2-point condensates outlined above can be generalized to (A^ > 2)-point con- 
densates by an auxiliary increase to A^ — 1 spacetime dimensions on which the functions / and g depend. 
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7.3.3 The gluonic field strength correlator and bilocal quark condensate on the lattice 

Here we review what is known from lattice simulations about the gauge invariant, gluonic field strength 
correlator which bears information on the nonlocal gluon condensate of Eq. (jll2|) . Appealing to 
Poincare, parity and time inversion invariance in Euclidean spacetime, the following parametrization 
was introduced in 



dDi{x'^) 

-\- ( X f^x i^5i,\ X ^X\5i^i^ + XyX\6^i^ ^u-^K^fiX ) Qx"^ ' ('^^^) 

To separate perturbative from nonperturbative contributions the scalar functions D and Di are usually 
[M MH M fitted as 
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V A3 



Di(x^) = Aiexp +-^exp . (123) 



The power-like behavior in Eqs. p23j) at small |a;| is believed to catch most of the perturbative physics, 
although it is known that partially summed perturbation theory may generate less divergent renormalon 
contributions at = as well, see Sec. (jS}- As we have seen, the coefficients of these contributions 
are ambiguous. We ignore this subtlety and assume that the purely exponential terms in Eqs. ()123p 
exclusively carry the nonperturbative information^^. We also note that due to the cusp of the exponential 
ansatz in Eq. ()123|1 . derivatives of the nonperturbative part have nonlogarithmic UV singularities at 
a; = that can only be defined in a nonperturbative UV regularization scheme. This is obviously in 
confiict with the well-tested idea that UV singularities, when treated in renormalization-group improved 
perturbation theory, lead to a logarithmic, anomalous scaling of averages over local operators. It does 
then make no sense to assume the exponential behavior of Eq. ()123jl down to arbitrarily small distances - 
after all the lattice resolution at which the field strength correlator is measured is finite, a^^ ~ 2 GeV. We 
will re-address this concern in Sees. 17. 3 ."il and rr3 . 71 We also would like to mention that in the framework 
of a nonrelativistic effective-field-theory formulation of QCD a decomposition of the correlator p22|) into 
color magnetic and color electric correlations yields an electric correlation length, which is smaller than 
the magnetic, one as a result of lattice measurements |252U250l I251j . See also |261j for a relation between 
the decay widths of heavy quarkonia on the one hand and electric and magnetic gluon correlators on 
the other hand. 

In a recent unquenched lattice simulation |246j (see also Ref. [247]), where the gluon field strength 
correlator was measured with a resolution of ~ (0.1 fm)~^ ^ 2 GeV between 3 and 8 lattice spacings, 
it was found that 

4^ ^ 9 and A-^ ^ 0.7 GeV. (124) 

Notice that the inverse correlation length is somewhat larger than the typical hadronization scale 
Aqcd- While the actual size of Ag and Ai depend quite strongly on whether quenched or unquenched 
simulations are carried out and also on the values for the light quark masses assumed, the ratio Ag/ Ai 
and the correlation length A^ were found to be quite stable |246j . The Di contribution can be neglected 
since Ag ^ Ai (see Eq. fll22|) ). and one may then write 

V'iA^''*!^;) = ^ ( Sf^JuX - S^x6^^ ) . g{\x\) (125) 
^^This is in contrast to |197l 11981 1199L EUni I2U5| where a Gaussian behavior was assumed. 
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This has the same tensor structure as the local condensate g^Jl^x^'^^ (0). It is in agreement with the local 
expansion of g^u^x which reads |il40j 

+ iS^p5au + S^pSaiy - - S^pS^a) , (126) 

where 




(127) 



and denotes a light flavor-singlet current. Using exact vacuum saturation for the four-quark oper- 
ator, asin = 0.7GeV) = 0.7, and the following instanton-calculus determined value of the condensate 
{9^ fabcGl^Gly^Gl^ = 0.045 GeV^ and (qq) = -(0.24 GeV)^, it can be shown that the tensor structure 

belonging to the local condensate dominates the remainder in {g'^G'^^j^dpdpG'^^j^j by a factor of about 
nine |7S] . 

In the case of the scalar, gauge invariant quark correlator the situation is similar |248j . For the 
longitudinal vector quark correlator defined as 

-^g(a:)7Mk,0]g (128) 

an about ten times smaller correlation length as compared to the scalar case was obtained in |248j with 
a four-flavor simulation (staggered fermions) at a lattice resolution of a^^ ~ 2 GeV. 



7.3.4 Model calculation of the nonperturbative shift in the ground-state energy of heavy 
quarkonium 

In this section we demonstrate in a model calculation for the nonperturbative shift of the ground-state 
[n'^^^^Lj = l^Si) energy E"'^ of the heavy quarkonia, where the "exact" result can be calculated, in de- 
pendence on the heavy-quark mass m how much better the delocalized operator expansion converges as 
compared to the local expansion [73j. Our investigation is not intended to represent a phenomenological 
study of non-perturbative effects in heavy quarkonium energy levels. 

Among the early applications of the OPE in QCD was the analysis of nonperturbative effects in heavy 
quarkonium systems |253| I254t I255j . Heavy quarkonium systems are nonrelativistic quark- ant iquark 
bound states for which the following hierarchy of the relevant physical scales m (heavy quark mass), 
mv (relative momentum), mw^ (kinetic energy) and Aqcd holds: 

m > mf > mv^ > Aqcd- (129) 

The spatial size of the quarkonium system ~ (mf )~^ is much smaller than the typical dynamical time 
scale ~ (mf ^)~^. In practice the last of the conditions in Eq. (jl29|) . which relates the vacuum correlation 
length with the quarkonium energy scale, is probably not satisfied for any known quarkonium state, 
not even for T mesons fl81l. Only for top-antitop quark threshold production condition ()129|) may be 
a viable assumption 2 57. . 

We adopt the local version of the multipole expansion (OPE) for the expansion in the ratios of 
the scales m, mv and mv'^. The resolution dependent expansion (DOE) is applied with respect to 
the ratio of the scales mv'^ and Aqcd- The former expansion amounts to the usual treatment of the 
dominant perturbative dynamics by means of a nonrelativistic two-body Schrodinger equation. The 
interaction with the nonperturbative vacuum is accounted for by two insertions of the local xE dipole 
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operator, E being the chromoelectric field. |253j The chain of VEV's of the two gluon operator with 
increasing numbers of covariant derivatives times powers of quark-antiquark octet propagators |253j . 
i.e. the expansion in A/mf^, is treated in the DOE. 

At leading order in the local multipole expansion with respect to the scales m, mv, and mv'^ the 
expression for the nonperturbative corrections to the ground state energy reads 



= / dt fit) git), (130) 
where 

m = ^ / ^ ^''"^''^ J J d'y H^) (xy) Go (x, y , - qo) Hv) , (i3i) 



with 

GoUy,--) = f:i2l + l)pJ^)GJx,y 

\ 771/ J ^ Q \ / \ 



xy / \ m 



Gi(x,y,--) = ^ i2kxy i2kyy e-'^^+y^ 

^3/2 



s=0 ys + l + l-^)is + 2l 



(pix) = —i=e 



-kx 



2 

k = -mas (132) 
3 

The term Go is the quark-antiquark octet Green-function |255j . and denotes the ground state wave 
function of the quarkonium system. The functions P„ and are Legendre and Laguerre polynomials, 
respectively. Since we neglect the spatial extension of the quarkonium system with respect to the 
interaction with the nonperturbative vacuum, the insertions of the xE operator probe only the temporal 
correlations in the vacuum. This effectively renders the problem one-dimensional. Notice that t is the 
Euclidean time. 

Since the spatial extension of the quarkonium system is neglected and the average time between 
interactions with the vacuum is of the order of the inverse kinetic energy, the characteristic width of 
the function / in Eq. ()13H) is of order (mf ^)~^ ~ (ma^)~^. The values of the first few local multipole 
moments /„(oo), which correspond to local Wilson coefficients, read 

/o(oo) = 1.65183!^, /2(oo) = 1.3130 f^. 



/4(oo) = 7.7570:^, /6(oo) = 130.492^, (133) 



, , 11 



/8(oo) = 4474.1 ^ , /lo(oo) = 262709.3 fp4 , • • • • 

The term /o(c)o) agrees with Ref. |254t I255j and /2(c)o) with Ref. |258j . The results for /n>2(oo) are new. 

We use a lattice-inspired model function git) for the nonperturbative gluonic field strength correlator 
of the form 



git) = 12Agexpi-Jt^ + Xl/Xg + l 



Ag = 0.04 GeV\ A^^ = 0.7GeV, (134) 



This model has an exponential large-time behavior according to Eq. p23|) and a smooth behavior for 
small t. The local dimension four gluon condensate in this model is 



- G%G;^) = ^ = 0.024 GeV^ . (135) 
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= 


= oo 


n = 


k'^/m 


m 
(GeV) 




k'^/m 
(MeV) 


(MeV) 


n 


fnQn 

(MeV) 


Ei=0 f'i9i 

(MeV) 


fnQn 

(MeV) 


J2i=0 fi9i 

(MeV) 


5 


0.39 


0.338 


24.8 





38.6 


38.6 


24.2 


24.2 










2 


-65.7 


-27.2 


-3.9 


20.3 










4 


832.7 


805.5 


12.1 


32.4 










6 


-35048.0 


-34242.4 


-43.1 


-10.8 


25 


0.23 


0.588 


12.6 





16.0 


16.0 


12.8 


12.8 










2 


-9.0 


7.0 


-1.2 


11.5 










4 


37.8 


44.8 


2.6 


14.1 










6 


-526.6 


-481.8 


-6.7 


7.4 


45 


0.19 


0.722 


4.9 





5.9 


5.9 


5.0 


5.0 










2 


-2.2 


3.7 


-0.4 


4.6 










4 


6.1 


9.8 


0.7 


5.3 










6 


-56.4 


-46.6 


-1.4 


3.8 


90 


0.17 


1.156 


1.05 





1.15 


1.15 


1.07 


1.07 










2 


-0.17 


0.98 


-0.04 


1.02 










4 


0.18 


1.16 


0.04 


1.07 










6 


-0.65 


0.51 


-0.07 


1.00 


175 


0.15 


1.750 


0.245 





0.258 


0.258 


0.249 


0.249 










2 


-0.016 


0.242 


-0.005 


0.244 










4 


0.008 


0.250 


0.003 


0.247 










6 


-0.012 


0.237 


-0.003 


0.244 



Table 1: Nonperturbative corrections Enp to the heavy quarkonium ground state level at leading order 
in the multipole expansion with respect to the scales m, mv and mv'^ for various quark masses m based 
on the model in Eq. ()134|) . Displayed are the exact result and the first few orders in the DOE for Q = oo 
and Q = k'^/m. The numbers are rounded off to units of 0.1, 0.01 or 0.001 MeV. 



The exact form of the model for the nonperturbative gluonic field strength correlator is not important 
for our purposes as long as the derivatives of g{t) at t = are well defined, see e.g. Refs. |18H 1259] for 
different model choices. Tab.^ shows the exact result in the model (jl34j) and the first four terms of 
the DOE of -E"^ for quark masses m = 5, 25, 45, 90, 175 GeV and for 1] = oo and Q = k'^/m. For each 
value of the quark mass the strong coupling has been fixed by the relation as = ctsik). Notice that the 
series all appear to be asymptotic, i.e. they are not convergent for any resolution. The local expansion 
{Q = oo) is badly behaved for small quark masses because for k'^/m < any local expansion is 
meaningless. In particular, for m = 5 GeV the subleading dimension-six term is already larger than the 
parametrically leading dimension-four term. For quark masses, where k'^/m > X~^, the local expansion 
is reasonably good. However, at the finite resolution Q = k'^/m, the size of higher order terms is 
considerably smaller than in the local expansion for all quark masses, and the series is apparently much 
better behaved. The size of the order-n term is suppressed by approximately a factor 2~" as compared 
to the order-n term in the local expansion. We see explicitly that terms in the series with larger n 
increase more quickly in magnitude in the local expansion scale as compared to expansion at finite 
resolution. One also observes that even in the case fc^/m < A~^, where the leading term of the local 
expansion overestimates the exact result, the leading term in the delocalized expansion for Q = k'^/m 
agrees with the exact result within a few percent. 
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For a realistic treatment of the nonperturbative contributions in the heavy quarkonium spectrum a 
model-independent analysis should be carried out. In addition, also higher orders in the local multipole 
expansion with respect to the ratios of scales m, mv and mv"^ should be taken into account, which have 
been neglected here. These corrections might be substantial, in particular for smaller quark masses. 

7.3.5 Running gluon condensate from the cc spectrum 

In this section we review an extraction of the running gluon condensate from the charmonium spectrum 
using moment sum rules [73J. We have already discussed this approach in Sec. ()4.2j) where an OPE was 
assumed (for an analysis involving operators at n = 8 see j72]. The nth moment Ain reads 

Ain = 7—5- — ^7— — — ^— J-. — 7— 1 1 + [pert, corrections] 



where 



47r2 (2n + 3)!! (4m2)" 



(n + 3)! 



+ 



(n-l)!(2n + 5) 9 (Ami) 



2 



5, 



(6) _ 2 {n + A)\ {3n^ + 8n - 5) 



45(n-l)!(2n + 5)(2n + 7) 9(4m2)3 
.(6) 8 (n + 2)!(n + 4)(3n3 + 47n2 + 2..,„ , . 



135 (r2-l)!(2n + 5)(2n + 7) 9 (4m2)3 

and {g'fG') = (g^r'^G^Gt^Gl^) = 0.045 GeV^ (instanton gas approximation |2T1), {g^f) = {g^J^) = 
-p4/3(47r)^Q;^(gg)^ (as(/i = 0.7 GeV) = 0.7 and (qq) = -(0.24 GeV)^, p = 1 ^ exact vacuum satura- 
tion), being the light flavor singlet current. We work with the ratio 

- (137) 



Mn- 



The extraction of the running of the gluon condensate is reliable^^ provided that subleading, dimension- 
six power corrections in Eq. ()136|) are much smaller than the dimension-six power correction 

-^6l^Hg'GD'G) (138) 

stemming from the local expansion of the dimension- four running gluon condensate. Table |21 shows that 
in the OPE for r„ this is indeed the case provided that n is sufficiently large (in practice n > 4). As 
expected. Table El indicates that for a value of the four-quark condensate twice the value obtained from 
exact vacuum saturation the convergence of the dimension-six part of the running gluon condensate 
towards that of the full OPE is slower than in the case of exact vacuum saturation. On the experimental 
side of the sum rules we have used the spectrum as compiled in jHl], see this reference for details. To 
compare the running gluon condensate as extracted from the data with the model expression (using 
the lattice fit of Eq. p23|) and the Ui = expression in Eq. pifi|) ) we have to make a choice for the 
resolution scale Q. The following physical arguments apply: For large n, i.e. in the nonrelativistic 
regime, the width of the short-distance function / is of the order of the quark cm. kinetic energy mw^. 



^•^Recall, that the gluon condensate perturbatively is a renormahzation-group invariant and thus docs not scale loga- 
rithmically. 
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n 


1 


2 


3 


4 


5 


6 


7 


8 


p=l 


0.13 


0.36 


0.56 


0.73 


0.87 


0.99 


1.08 


1.16 


p = 2 


0.09 


0.25 


0.42 


0.57 


0.69 


0.80 


0.89 


0.97 



Table 2: Ratio of the local dimension-six contributions contained in the running gluon condensate and 
in the full OPE for r„ as a function of n. 




Figure 25: The running gluon condensate as a function of n when extracted from the ratio 
of moments r„ for MS charm-quark masses rfic{jnc) = 1-23 (white triangles), 1.24 (black 
stars), 1.25 (white squares) and 1.26 GeV (black triangles). The area between the upper 
and lower symbols represents the uncertainties. The thick solid line indicates the running 
gluon condensate as it is obtained from the lattice-fitted ansatz in Eqs. ()123p . Plot taken 
from 



which scales like m/n because the average quark velocity in the n-th moment scales like ^l\fn |260j . 
For small n, on the other hand, the relevant short-distance scale is just the quark mass. We take this 
as a guideline to use f2 = Imd n as the relevant resolution scale in the model expression 




(139) 



for the running gluon condensate (A~^ = 0.7 GeV and Ag = 0.04 which corresponds to {{as/ 7r)G'^)iat(c)o) = 
0.024 GeV^). In Fig. ()25j) the results of the data extraction and the model calculation are shown. The 
n-dependence of the lattice-inspired running gluon condensate^"^ and the result obtained from the char- 
monium moment sum rules are consistent for larger n. For small n no conclusive statement can be 

^^Here and in Sec. l7.3.5l a steepening of the model curve (a stronger running) at large momenta (small n) will occur if 
the lattice fit of Eq. (|123|l is considered as the result of a finite-resolution measurement, see Sec. l7.3.7l 
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made, recall Table|21and the fact that for small n the sensitivity to the error in the continuum region 
of the spectrum is enhanced. 



7.3.6 Running gluon condensate from the r decay spectrum 

In this section we review an independent extraction of the running gluon condensate from the spectral 
function of the V + A channel |^ as it was measured in r decays by Aleph [179j and Opal |18Uj at LEP. 
In this channel the OPE of the associated current-current correlator i J dlce*^^ ^Tj^(x)jj^(0)^ (with 

currents j^^^ = M7/i(l ± 75)*^) is dominated by the gluon condensate |176[ 11771 1178j . the dimension- 
six power corrections that are not due to the local expansion of the running gluon condensate are 
suppressed. We use a sum rule for the cutoff independent Adler function 



Dm^~Q' '''\^r =- '7 - (140) 



dl V^+^{Q^) _ , Imnv+^(g) 

W ~ ~ Jo (s + g2)2 

For the V + A spectral function we have used the Aleph measurement [179] in the resonance region up 
to 2.2 GeV^. For the continuum region above 2.2 GeV^ 3-loop perturbation theory, which we also used 
on the OPE side, was employed (with as{Mz) = 0.118), and we have set the renormalization scale /i 
equal to Q. The Wilson coefficient for (^^G^^ was taken into account up to order jl31j . It can be 
read off from 

w^^(i-^^)e«v-- 

With the same values for the local condensates as in Sec. 17. 3. 51 the evaluation of the dimension-six 
contribution in the OPE and the dimension-six contribution of the local expansion of the running 
gluon condensate yields a comparable magnitude and equal signs which gave us a phenomenological 
justification for the extraction of the running gluon condensate in the V + A channel. Fig. 1^ shows the 
result. The uncertainties are due to the experimental errors in the spectral function and a variation 
of the renormalization scale fi in the range Q ± 0.25 GeV. The analysis is restricted to the range 
1 GeV < Q < 2 GeV because for Q < 1 GeV perturbation theory becomes unreliable and for Q > 2 GeV 
the experimentally unknown part of the spectral function at s > 2.2 GeV^ is being probed. The thick 
black line in Fig.l^ shows the lattice-inspired model for the running gluon condensate of Eq. ()139|) for 
Q = Q. It is consistent with the phenomenological extraction and an increasing function of Q. Since the 
Q dependence of the model is rather weak for 1 GeV < Q <2 GeV and the error band of the extraction 
rather large it is not possible to draw a more quantitative conclusion at present. 



7.3.7 Euclidean position-space V±A correlators at short distance 

In Sec. 17.3. !2] a delocalized version of the OPE was obtained by projections of a perturbatively calcu- 
lable short-distance function f{x) and a nonperturbative, long-distance function g{x'^) on a resolution 
dependently "rotated" basis in dual space, see Eq. (jll4|) . Thereby, the function gi^x"^) is assumed to be 
determined at infinite resolution. However, the only so-far available first-principle approach to 5'(x^) is 
a lattice calculation which can only be performed at a finite resolution, Q ~ a~^. 

An alternative approach to running condensates than the one in Sec. 17. 3. 31 is to view the OPE as 
the usual local expansion but now involving averages over local composite operators at a resolution 
Q Q. These operator averages are, besides their usual perturbative evolution, nonperturbatively 
coarse-grained in a gauge invariant way |145j . In the case of a composite with two fundamental operators 
like the quark or the gluon condensate nonperturbative coarse graining is performed by integrating out 
length scales between and (f2 — dQ)~^, which are associated with a loss in resolution of Q'^/dQ, in 
an fl dependent, gauge invariant and nonperturbative correlation function g{x, fl), using a sharply cut 
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Figure 26: The running gluon condensate as a function of Q when extracted from the Adler 
function. The grey area represents the allowed region using perturbation theory at 0{al) 
and the striped region using perturbation theory at 0{al). The thick sohd hne denotes the 
running gluon condensate as obtained from the lattice-fitted ansatz in Eqs. (|123|) for Q = Q. 
Plot taken from [7H] . 



off spherical well 

^^eii/dn-\x\) (142) 

as a weight function. If we assume self-similarity of g{x, Q), that is, an exponential form 

gix,n)=Ain)exp[-\x\/Xg], (143) 

where only the coefficient^^ A{Q) depends on Q, then the following evolution equation for A{Q) is easily 
derived |145j 

A^(n) = ±n-M(n). (144) 

The solution of Eq. ()144j) corresponding to the running gluon condensate would then take the form 

1 1 



;^Gn(^]) = '4 

TV I TT^ 



(145) 



when the lattice fit of Eq. ()123|) is used as a model with fliat = = 2 GeV. 

Following |147lj we now use dependences like Eq. ()142|) to investigate the effect of running condensates 
on the |x| dependence of the Euclidean V ± A position-space correlators when they are expanded into 
a DOE in the chiral limit [TTT] 



3 



\X\ 



^^Interpreting Ag ^ as the mass of the lowest intermediate hadronic state reached by the gauge invariant correlation, 
see for example jMS) , the correlation length Xg is viewed as an observable which does not depend on the resolution. 
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X [GeV^'l 



X [GeV-'l 



Fi gure 27: Ry^A cis a function of distance x in the chiral limit. The sohd hne corresponds 
to the DOE with an assumed factorization of four-point quark correlators at -D = 6 into 
two-point, scalar quark correlators, the dashed line to the OPE. Crosses depict the result of 
the instanton liquid calculation of ^126^ which is taken from |121j . Plot taken from |147j . 



Rv+a{\x\ ~ 




6 



(146) 



in the V ± A channel. There is no perturbative and gluon-condensate contribution in Rv-a, the 
correlator is extremely sensitive to chiral symmetry breaking. For the treatment of the four-point quark 
correlators at D = 6 we assume a factorization into scalar two-point quark correlators. In analogy to 
the case of the field strength correlator a lattice measurement of the scalar two-point quark correlator 
j248j with Np = 4 staggered fermions and a quark mass am = 0.01 at ~ 2 GeV yields 



Substituting Q = l/\x\ in Eq. ()146|) and working with the parameters of Eqs. p47p . p34|) . and a fixed 
value^^ of as{Qiat) = 0.2 |249j . an |x| dependence as depicted in Fig. |2Z|is obtained. In addition to 
a much better agreement (in comparison to the OPE) with the result obtained in the instanton liquid 
|126j the result of the DOE calculation almost perfectly agrees with the result of a quenched lattice 
calculation obtained using an overlap action |121j . The result of an extraction of -Rvi^d^^l) from the 
r-decay data, see for example |126j . slightly overshoots the DOE and lattice results which were obtained 
in the chiral limit see Fig.OHl 

7.3.8 Mesonic spectra from the DOE? 

The good agreement that we have found in Sec. 17. 3. 71 between Euclidean V ± A correlators when ex- 
panded into an OPE with running condensates and when computed on a lattice lead us to investigate to 
what extent one can predict properties of the spectral functions of light-quark channels within the reso- 
nance region from the (practical) OPE with running condensates by analytical continuation to time-like 
momenta, = — (s -|- ie) or Q = —i-\/s, (s > 0), and by taking the imaginary part afterwards |146j . 
For a contribution of the form 



Xq = 3.1 GeV 



1 



Aq{a 



Vtiat ~ 2 GeV) = (0.212 GeV)^ . 



(147) 




(148) 



The running of ag almost cancels the log-powers at D = 6 in Eq. (|146|l . 
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Figure 28: Rv±A as a function of distance x. The dashed 
hues indicate the errors in the experimental determination 
of the respective spectral functions by Aleph. The squares 
show the results obtained in the model of a random instan- 
ton liquid. The diamonds are the OPE prediction using 
realistic light-quark masses. Plot taken from jl2bj . 



the associated contribution to the spectral function reads 



lat 



exp 



5A4 Vtiat 



sm 



(149) 



Notice that oscillations only start if \fs becomes smaller than the inverse correlation length Aj^. As 
in Sec. 17. 3. 71 a factorization of four-quark correlators into scalar two-quark correlators within the usual 
vacuum saturation hypothesis of the local case is assumed. All running condensates thus have the form 
as in Eq. ()149|) . the effective correlation length for running condensates of D = 6 is Aq/2. The p, ai, tt, 
and meson channels where investigated in |146j . and the following two sets of parameters were used: 



SetA : 
Sets : 



3.1 GeV"^ 
1.7 GeV"^ 



2GeV) 
2GeV) 



(0.212 GeV)^ 
0.015 {GeVy 
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(iiB) (iiiB) (ivB) 



Figure 29: The spectral functions Imn(Q^ = —s — ie), (s > 0) for sets (A) and (B), where 
(i), (u), (iii), and (iv) correspond to the p, ai, vr, and channels, respectively. The unit of 
s is (GeV)2. Taken from [HH] . 



Ag = 0.3 GeV"^ , Ag{a-^ = Qiat ~ 2 GeV) = (0.212 GeV)^ , 

Xg = 1.7 GeV"^ , Ag{a-^ = ^i^t ~ 2 GeV) = 0.015 (GeV)^ . (150) 

The quark correlation length Xg in Set A corresponds to that of the scalar quark correlator, the one 
in Set B to the longitudinal vector quark correlator as obtained in |248j . Judging from Fig.EHl the 
result obtained with the small fermionic correlation length in Set B is more realistic but still far off the 
experimentally observed behavior. The most disturbing feature of the spectra in Fig. 1211 is the fact that 
they become negative at low values of s. For a discussion of the moments of the spectra in Fig.l^see 
|146j . This is unacceptable and certainly has to do with the incompleteness and hence poor practical 
convergence of the expansion at low momenta. 



7.3.9 Speculations on convergence properties of expansions with nonperturbatively run- 
ning operators 

How does the partial resummation of operators involving powers of covariant derivatives, which leads 
to the nonperturbative running of the condensates in the usual OPE, possibly affect the "convergence" 
properties of the modified operator expansion? Let us give some (admittedly speculative) arguments. 

First of all, one notices the difference between the exponential(-like) dependences of Eq. (|145|) (or 
Eq. ()139|) ) with the situation in the instanton model in Sec. 16. 21 where exponentially small terms of the 
form exp[ ^^^ ] were observed. As it was argued in |142j the latter exponentials may occur in the 

exact result to cure the ambiguities of the 1/Q^ expansion which may arise due to the factorially-in- 
D rising coefficients where D roughly refers to the power in 1/Q. A possibility is that this factorial 
divergence of coefficients would already be present in a modification OPE of the OPE where operators 
containing powers of covariant derivatives are omitted^^ - the factorially-in-D rising coefficients would 
then simply arise from the combinatorial variety of composites involving powers of quark and gluon 
field operators. So already the OPE and, more generally, each subseries of the OPE with a given, fixed 

^^We do not want to apply the equations of motion and Bianchi identities which reduce these to operators without 
powers of covariant derivatives. 
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power of covariant derivatives would be an asymptotic expansion which could be made unambiguous 
by adding exponentially small terms of the form exp[j^ — l 



^QCD' 



Second, in the OPE we may resum a part of the powers-of-covariant-derivative series associated with 
each operator in OPE by applying our methods of Sees. 17.3.'^ or I7.3.'7I For an operator of dimension D 
the associated contribution after partial resummation is roughly of the form 



lat 



exp 




lat , 



(151) 



where A„ denotes an effective correlation length. The expression in Eq. ljlSip has a maximum at Q 
inXr,)^^ with value 



Pn = (An A„ exp[-l] n)"' exp 



n 



lat 



(152) 



where A„ = (An{^iat)Y^^'^^ ■ As can be motivated from the example of factorizing the four-quark 
correlations into two-quark correlations in Sec. 17.3.71 it is likely that A„ is a decreasing function of n. If 
we assume a fall-off as A„ = \/n^^~^\ (1 > e > 0) then the position of the maximum of the expression 



in Eq. p51|) decreases as Qr 



X- 



n 



The value of the maximum in this case is 



Pn = A„An" exp[-l + 



n 



xn 



lat 



(153) 



Assuming that A„ does not depend on n, A„ = Aqcd, as suggested by the phenomenology of practical 
OPEs, the critical mass dimension Uc from which on maxima explode is given in implicit form as 



AqcdXuI exp[-l + 



xn 



lat 



(154) 



A truncation of the OPE with running condensates at some n < ric would mean that physics below 
the maximum Q„ can not be described, a truncation at n > ric does not make sense since due to the 
rapid increase of maxima aX n > Uc the expansion does not approximate anymore. Provided our above 
assumptions are met, another source of asymptotic behavior is identified. Let us give some numerical 
examples. For e = 0.5 and Aqcd = 0.4 GeV (the rounded) value of was calculated as a function of 
A in [Hg as 



A = 3GeV"^ ^ n, 
X = IGeV-^ rir 



4 ; A = 2 GeV-^ 
40; A = 0.5GeV-^ 



rir 



10; 
158, 



In conclusion, we have argued that exponential-like behavior of the form 



exp I 



Q 



A, 



and 



exp 



QCD 



-A 



-1 



n 



lat . 



(155) 



(156) 



may coexist in an improved version of the OPE. The former may assure the uniqueness of the asymptotic 
1/Q expansion involving operators with a given, fixed number of covariant derivatives |142j while the 
latter arises from a partial summation of the expansion in powers of covariant derivatives acting on a 
given, fixed composition of powers of quark and gluon operators. 



8 Summary 

In this article we have presented a review of theoretical approaches to the indentification of the mech- 
anisms leading to the violation of local quark-hadron duality in QCD. We have started our discussion 
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with a mini-review on QCD sum rules in (axial)vector-meson channels in vacuum, at finite temperature, 
and at finite baryon density. The status of renormalon singularities as possible sources for power correc- 
tions in the OPE was very bricfiy addressed. We then have reviewed two theoretical model approaches 
to the calculation of current correlators: the instanton-gas model and the 't Hooft model. The largest 
part of the review was devoted to a discussion of phenomenlogical evidence for the violation of quark- 
hadron duality. We have gathered indications that nonperturbatively nonlocal effects are insufficiently 
incorporated in a truncated local expansion if nonlocal quantities like parton distribution amplitudes, 
meson form factors, and hadron transition amplitudes are to be predicted. Phenomenologically moti- 
vated solutions to this shortcoming were discussed. Possible imphcations for OPE-based predictions of 
hadron spectra were speculated upon. 
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